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ABSTRACT 

This  work  investigates  the  fundamental  relationships  of  wave  loading  on  cross  deck 
structures  for  trimaran  vessels.  In  contrast  with  a  monohull  ship,  trimaran  vessels  experience 
several  possible  structural  loading  cases  including:  longitudinal  bending,  transverse  bending, 
torsional  bending,  spreading  and  squeezing  of  hulls,  inner  and  outer  hull  slam  pressures,  wet 
deck  slam  pressures,  loading  from  ship’s  motions,  and  whipping  of  slender  hulls.  This  work 
investigates  wave  loading  cases  that  result  in  transverse  and  torsional  bending  of  the  cross  deck 
structure. 

The  wave  loading  cases  investigated  include:  side  hull  troughing  and  cresting  in 
longitudinal  waves;  side  hull  torsion  in  longitudinal  waves,  and  transverse  hogging  and  sagging. 
For  each  of  these  load  cases,  a  design  load  using  a  fully  statistical  sea  state  was  derived  using  an 
analytical  model  of  a  trimaran  represented  by  rigidly  connected  box  barges.  The  design  loadings 
with  a  reliability  index  of  5  for  almost  500  trimaran  configurations  were  calculated  varying  main 
hull  length,  side  hull  length,  side  hull  transverse  placement,  and  side  hull  longitudinal  placement. 
The  design  loadings  were  curve  fit  to  a  fourth  order  polynomial  in  the  three  independent 
variables. 

The  load  predictions  of  the  analytical  box  model  of  a  trimaran  were  applied  to  a  trimaran 
vessel  with  a  realistic  hull  form  using  the  finite  element  ship  structural  analysis  program 
MAESTRO.  Given  the  number  of  approximations  and  assumptions  in  the  analytical  model,  the 
forces  predicted  by  analytical  model  agreed  closely  with  the  finite  element  model’s  results. 

The  fitted  curve  of  design  loadings  allows  an  initial  design  stage  loading  estimate  for 
cross  deck  structural  loading,  given  general  characteristics  of  length  and  spacing  of  a  trimaran’s 
hulls.  This  estimate  of  structural  loading  combined  with  other  characteristics  of  good  trimaran 
design  including  stability,  roll,  and  resistance  characteristics  will  aid  in  optimizing  an  overall 
trimaran  ship  design. 
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Chapter  1  Introduction 


Over  recent  decades  a  growing  demand  for  higher  speed  ships  has  led  to  the  development 
of  several  new  hull  form  concepts.  These  hull  forms  include  catamarans,  surface  effects  (SES) 
ships,  small  waterplane  twin-hull  (SWATH)  ships,  pentamaran,  and  trimaran  ships  to  name  a 
few.  Of  these  ‘new’  hull  forms  the  multi-hull  vessel’s  origins  can  be  traced  to  back  several 
centuries  to  many  seagoing  peoples  with  outrigger  canoes.  In  more  recent  times  the  multi-hull 
ships  have  been  used  in  several  racing,  pleasure,  and  commercial  vessels. 

The  benefits  of  the  trimaran  hull  form  have  been  studied  extensively  over  the  past  several 
years  at  the  University  College  London  (UCL)  [1].  In  2000  the  RV  TRITON,  Figure  1  below,  a 
trimaran  demonstrator  project  for  the  United  Kingdom’s  Royal  Navy,  was  launched  to  test  the 
trimaran  hull  form  [2]. 


Figure  1  -  RY  Triton  [2] 

While  much  work  in  the  area  of  trimaran  hull  form  design  has  been  accomplished  by  the 
UCL  studies,  the  structural  loading  experienced  by  the  cross-deck  structure  of  the  trimaran  hull 
form  is  still  largely  unknown.  Classification  societies  such  as  the  American  Bureau  of  Shipping 
(ABS)  [3]  or  Det  Norske  Veritas  (DNV)  [4]  currently  have  design  codes  for  a  traditional 
monohull’s  design  loadings  in  terms  of  the  ship’s  relevant  dimensions  and  a  sea-state  coefficient 
dependent  on  ships  dimensions  to  provide  sufficient  design  margin  for  the  life  of  a  ship. 
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However,  currently  no  such  design  codes  exists  for  trimarans  in  general  or  specifically  their 
cross-deck  structures.  The  goal  of  this  work  is  to  be  able  to  state  similar  design  loadings  based 
on  the  relevant  parameters  of  a  trimaran  ship. 

The  trimaran’s  structural  loading  will  depend  strongly  on  the  longitudinal  extent  as  well  as 
the  longitudinal  and  transverse  location  of  the  outer  hulls.  These  parameters  of  the  outer  hulls, 
by  necessity,  depend  on  the  operational  requirements  and  uses  of  the  ship.  In  contrast  with  a 
conventional  mono-hull  ship,  the  structural  loading  of  the  trimaran  potentially  involves  several 
additional  loading  cases  not  experienced  by  mono-hull  ships.  These  loading  cases  may  include 
longitudinal  bending,  transverse  bending,  torsional  bending,  spreading  and  squeezing  of  hulls, 
inner  and  outer  hull  slam  pressures,  wet  deck  slam  pressures,  loading  from  ship’s  motions,  and 
whipping  of  slender  hulls. 

1.1  Thesis  and  general  approach 

Of  the  above-mentioned  loading  cases  this  work  will  focus  on  the  structural  loading  in  the 
longitudinal,  transverse,  and  torsional  wave  loading  cases  that  affect  the  cross-structure  between 
the  main  and  outer  hulls  of  a  trimaran.  As  previously  mentioned,  the  loading  of  the  cross 
structure  between  the  hulls  of  a  trimaran  will  depend  strongly  on  the  longitudinal  extent  and 
location  of  the  trimaran’s  side  hulls.  This  work  will  quantify  the  effects  of  various  outer  hull 
placements  and  sizes  on  the  trimaran’s  cross-deck  structural  loading.  Once  the  structural  loading 
as  functions  of  placement  and  size  are  determined,  this  information  could  be  used  in  conjunction 
with  other  trimaran  design  parameters  such  as  stability,  roll,  and  resistance  characteristics  to 
optimize  an  overall  trimaran  ship  design. 

The  approach  for  this  study  will  be  to  first  determine  analytical  approximations  to  the 
trimaran  cross-deck  structural  loading  using  simple  and  symmetric  box -type  hull  shapes  for  each 
of  the  three  hulls.  These  analytical  solutions  will  account  for  worst-case  statistical  sea  state  for 
various  outer  hull  placement  and  size.  Next  these  analytical  approximations  will  be  compared 
with  quasi-finite  element  solutions  obtained  using  the  ship  structural  design  program  Method  for 
Analysis  Evaluation  and  Structural  Optimization  (MAESTRO)  on  hull  forms  more  closely 
resembling  actual  hulls  of  trimarans.  Comparison  of  these  two  approaches  will  provide  an 
estimate  the  validity  of  the  analytical  approximations  for  various  combinations  of  sea  state  and 
different  trimaran  configurations  of  placement  and  size  of  the  outer  hulls. 


1.2  Background 


The  openly  available  previous  work  in  the  area  of  trimaran  structural  design  to  date  has 
been  fairly  limited  in  scope.  In  general,  the  previous  studies  have  assumed  a  worst  case 
deterministic  loading  level  in  a  few  loading  cases  for  specific  designs  of  trimarans.  The  foci  of 
these  studies  have  been  to  investigate  the  cross-structure  contribution  to  primary  hull  bending 
and  transverse  cross-structure  bending  in  rolling  conditions.  The  results  of  these  studies  are 
included  in  references  [5]  through  [10]. 

The  first  of  these  studies  examined  the  contribution  from  the  side  hull  and  cross  deck 
structure  to  overall  ship  structural  performance  in  longitudinal  bending  for  both  hogging  and 
sagging  [5].  The  next  study  involving  trimaran  structures  focused  on  performing  a  detailed 
structural  point  design  including  scantling  sizes  to  estimate  the  structural  weight  fraction  of  a 
specific  frigate-sized  trimaran  as  compared  to  a  monohull  [6].  The  third  study  again  investigated 
the  contribution  of  side  hull  and  cross  deck  structure  to  the  resistance  of  primary  hull  bending  in 
hogging  and  sagging  [7].  The  next  study  [8]  involved  a  reevaluation  of  the  contribution  to 
primary  bending  resistance  from  the  side  hulls  and  cross  structure  contained  in  references  [5] 
and  [7]  with  a  more  refined  model  of  the  trimaran  hull  form. 

While  references  [5]  through  [8]  dealt  primarily  with  longitudinal  loading  in  primary 
bending,  reference  [9]  investigated  the  trimaran’s  structural  behavior  under  transverse  loading 
due  to  rolling  the  trimaran’s  side  hulls  to  complete  submergence  or  broach.  The  last  study  in 
trimaran  structural  design  [10]  investigates  low  weight  alternatives  to  cross  structure  in  loading 
condition  cause  by  a  ship  rolling  its  hulls  to  complete  submergence  or  broach. 

The  previous  work  on  trimaran  structural  design  to  this  point  has  only  scratched  the 
surface  of  the  possible  relevant  structural  design  issues  of  this  new  hull  form.  While  the 
previous  work  has  sought  to  characterize  a  few  deterministic  loading  scenarios  for  a  small  range 
of  ships,  this  work  will  characterize  the  structural  wave  loading  of  the  trimaran’s  cross-deck 
structure  subjected  to  a  statistical  sea  state  for  a  variety  of  trimaran  configurations  and  sizes. 


Chapter  2  General  Trimaran  Characteristics 

The  trimaran  hull  form  has  been  studied  extensively  over  the  past  several  years  at  the 
University  College  London.  The  results  of  these  studies  have  produced  several  variants  of  the 
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trimaran  ship.  Some  of  the  relevant  characteristics  of  those  designs  are  shown  below  in  Table  1. 


Table  1  -  Principal  characteristics  of  UCL  trimaran  studies  [1] 


Offshore 

Patrol 

Vessel 

1 

mi 

HzSSfl 

Corvette 

AAW 

Destroyer 

Cruise 

Liner 

Small 

Carrier 

Dispfacement  (tonne) 

234 

514 

1117.6 

ITBl 

1350 

1777 

4200 

4300 

4978 

9050 

11850 

16657 

61.04 

78.8 

98 

105 

120 

112 

154.7 

156.8 

168.6 

192 

191.5 

231.6 

10.85 

13.7 

22.5 

19.2 

25 

20 

27.5 

25.9 

25 

28 

ES 

43 

4.3 

8.5 

9 

8.5 

8 

8.85 

10.23 

12.1 

11.1 

13.2 

23.35 

23.5 

Uair  Hull  LWL  (m) 

59.8 

76.8 

91 

99 

115 

106.7 

148.7 

149.8 

151.3 

178.3 

177.2 

220 

Main  Hull  Beam  WL  (m) 

4.2 

4.2 

6.848 

6.8 

6.5 

8.5 

10.4 

10.8 

10.8 

13 

13.5 

14.5 

Main  Hull  Draft  WL  (m) 

2.1 

3.4 

3.2 

3.4 

3.2 

4.25 

5.2 

5.3 

4.8 

6.4 

8.74 

8 

Side  Hull  Displacement 

4.20% 

3.10% 

3.70% 

4.00% 

3.80% 

4.30% 

5.50% 

3.70% 

4.70% 

3.00% 

5.00% 

6.80% 

Side  Hull  LWL  (m) 

19.9 

28 

34.2 

35 

30 

50 

36 

56.9 

65 

71.3 

65.2 

82 

Side  Hull  Beam  WL  (m) 

1.06 

0.74 

1.45 

1.5 

2 

2.7 

3 

2 

2.5 

2.8 

3.65 

4 

Side  Hull  Draft  (m) 

0.9 

2.1 

2.31 

2 

1.5 

1.35 

3.6 

2.8 

2.7 

2.6 

4.37 

6.5 

Max  Speed  (knots) 

25 

25 

20 

38 

36 

30 

28 

28 

28 

26 

18 

27 

Ps  (MW) 

2.14 

4.3 

4 

20 

20 

20 

24 

26 

29 

31.5 

16.8 

70 

The  general  form  of  the  trimaran  in  these  balanced  trimaran  ship  designs  is  determined  by 
several  factors  including  stability,  roll,  and  resistance  characteristics  of  a  trimaran  ship.  In 
general  the  length  of  the  side  hulls  is  determined  by  intact  and  damaged  stability  requirements 
[1],  The  motion  characteristics  for  trimarans  have  been  examined  [1].  Another  study  indicates 
that  outer  hull  waterplane  area  strongly  affects  roll  motion  [1 1].  Favorable  resistance 
characteristics  of  trimaran  outer  hull  placement  have  also  been  studied  [12]. 


From  the  trimaran  characteristics  of  the  UCL  designs  in  Table  1,  other  important 
characteristics  of  balanced  trimaran  ship  designs  were  calculated  or  derived  with  the  results 
shown  below  in  Table  2. 


Table  2  -  Additional  characteristics  of  UCL  trimaran  studies 


Small 

Support 

Vessel 

Offshore 

Patrol 

Vessel 

RV 

TRITON 

Fast 

Ferry 

Canadian 

Ferry 

Corvette 

ASW 

Frigate 

ASW 

Frigate 

AAW 

Destroyer 

Cruise 

Liner 

LPH 

Small 

Carrier 

Main  Hull  Relations 

Main  Hull  UB 

14.24 

18.29 

13.29 

14.56 

17.69 

12.55 

14.30 

13.87 

14.01 

13.72 

13.13 

15.17 

Main  Hull  B/T 

2.00 

1.24 

2.14 

2.00 

2.03 

2.00 

2.00 

2.04 

2.25 

2.03 

1.54 

1.81 

Main  Hull  D/T 

2.05 

|  ■ 

2.81 

2.50 

2.50 

2.08 

1.97 

2.28 

2.31 

2.06 

2.67 

2.94 

Main-Side  Hull  Relations 

Side/Main  Hull  L/L 

0.33 

0.36 

0.38 

0.35 

0.26 

0.47 

0.24 

0.38 

0.43 

0.40 

0.37 

0.37 

Side/Main  Hull  T/T 

0.43 

0.62 

0.72 

0.59 

0.47 

0.32 

0.69 

0.53 

0.56 

0.41 

0.50 

0.81 

Side/Main  Hull  B/B 

0.25 

0.18 

0.21 

0.31 

0.32 

0.29 

0.19 

0.23 

0.22 

0.27 

0.28 

1.17 

1.54 

1.54 

1.30 

1.77 

1.02 

1.18 

1.11 

1.04 

0.97 

1.35 

1.34 

Side  Hull  Relations 

Side  Hull  L/B 

18.77 

37.84 

23.59 

23.33 

15.00 

18.52 

12.00 

28.45 

26.00 

25.46 

17.86 

20.50 

Side  Hull  B/T 

1.18 

0.35 

0.63 

0.75 

1.33 

2.00 

0.83 

0.71 

0.93 

1.08 

0.84 

0.62 

The  trimaran  characteristics  shown  in  Table  2  will  be  used  later  in  section  4.1  for 
development  of  the  analytical  trimaran  model  to  set  proper  limits  governing  the  range  of 
characteristics  for  investigation  of  the  applicable  trimaran  design  space. 


Chapter  3  Theory 

The  most  accurate  way  to  determine  a  ship’s  structural  loading  is  to  perform  a  fully 
dynamic  analysis  of  the  ship  in  a  completely  statistical  sea  state  accounting  for  added  mass  and 
damping  of  the  ocean  in  relation  with  the  ship’s  motions  with  provisions  to  add  the  effects  of  the 
ship’s  forward  speed  and  heading.  However,  even  with  today’s  advanced  computational 
capabilities,  this  level  of  analysis  is  prohibitive  for  an  investigation  of  the  basic  structural  loading 
attributes  of  a  new  hull  form  such  as  the  trimaran.  Therefore,  the  emphasis  of  this  work  will 
study  a  static  analysis  of  the  wave  loading  of  a  trimaran,  which  often  is  sufficient  as  a  first 
estimate  of  structural  loading  [13]. 

The  cross-deck  structural  wave  loading  of  the  trimaran  is  affected  by  two  major 
considerations:  the  longitudinal  and  transverse  placement  of  the  outer  hulls,  and  the  waterplane 
area  (length  and  beam)  of  the  outer  hulls.  How  each  of  these  outer  hull  characteristics  affect 
each  loading  case  is  detailed  below.  ■ 

3.1  Wave  Loading  Cases 

;  '  __  ■  :  *■  ,  ’  ’  . 

3.1.1  Longitudinal  Bending 


m 


The  longitudinal  bending  loading  of  a  traditional  monohull  ship  is  generally  characterized 
by  two  loading  cases:  hogging  and  sagging.  Similarly,  the  longitudinal  bending  loads  exerted 
on  the  trimaran’s  cross-deck  structure  might  be  expected  to  be  described  in  terms  of  the  hogging 
and  sagging  cases.  However,  in  the  case  of  a  trimaran  there  can  be  two  cases  of  hogging  or 
sagging.  These  cases  include  hogging  or  sagging  of  the  main  hull  as  well  as  that  of  the  outer 
hulls.  Longitudinal  bending  loads  in  hogging  and  sagging  of  the  side  hulls  was  investigated  in 
[8]  and  found  to  be  virtually  insignificant.  While  hogging  and  sagging  of  the  main  hull  must  be 
designed  for  in  the  overall  structural  design  of  a  trimaran,  the  primary  bending  conditions  of  the 
main  hull  of  a  trimaran  is  not  significantly  different  than  that  of  a  conventional  monohull  ship. 
Additionally,  the  contribution  to  resistance  of  primary  bending  of  the  ship  due  to  the  addition  of 
the  trimaran’s  cross-structure  was  investigated  with  the  results  in  references  [5]  through  [8],  and 
hence  the  main  hull  longitudinal  hogging  and  sagging  problem  is  not  further  investigated  in 
depth  in  this  work. 


# 


# 


# 
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3.1.2  Coincident  Side  Hull  Troughing  and  Cresting  in  Longitudinal  Waves 

While  the  longitudinal  bending  forces  incurred  from  the  side  hulls  experiencing  a  wave 
of  a  wavelength  that  can  produce  hogging  or  sagging  on  the  side  hulls  are  not  significant  [8],  the 
possibility  of  the  wave  crest  or  trough  coinciding  with  the  side  hull  length  can  create  a  significant 
vertical  structural  loading.  This  loading  is  due  to  the  coincidence  of  location  of  the  side  hulls 
with  a  trough  or  a  crest  of  a  wave.  While  specific  wavelengths  can  create  a  significant  vertical 
force  on  the  side  hulls  due  to  broaching  or  immersing  the  side  hulls,  the  overall  structural 


response  of  the  entire  ship  is  minimal  compared  with  the  wavelengths  associated  with  primary 
bending  in  the  trimaran’s  main  hull  hogging  or  sagging.  These  conditions  are  shown  below  in 
Figure  2  through  Figure  5  with  the  hulls  of  the  trimaran  represented  by  box  barges. 


Figure  2  -  Longitudinal  Cresting  with  Side  Hulls  Aft 


Figure  2  through  Figure  5  show  profile  views  of  various  possible  configurations  of 
coincident  side  hull  troughing  and  cresting  for  a  trimaran  of  typical  proportions  as  determined 
from  Table  2.  The  dotted  lines  show  the  still  water  position  of  the  trimaran  and  waterline,  the 
dashed  line  shows  the  wave,  and  the  solid  lines  show  the  trimaran’s  response  to  the  wave. 


Figure  3  -  Longitudinal  Cresting  with  Side  Hulls  Amidships 


Figure  4  -  Longitudinal  Troughing  with  Side  Hulls  Aft 


Figure  5  -  Longitudinal  Troughing  with  Side  Hulls  Amidships 

In  contrast  with  the  traditional  primary  bending  of  hogging  and  sagging,  coincident  side 
hull  troughing  and  cresting  creates  a  vertical  loading  that  increases  and  decreases  the 
contribution  of  side  hull  buoyancy.  For  example,  as  can  be  seen  in  Figure  5  above,  the  main  hull 
experiences  a  partial  sagging  condition,  while  the  trough  of  the  wave  drops  out  from  beneath  the 
side  hulls  creating  a  vertical  downward  force  on  the  cross-structure  due  to  the  loss  of  buoyancy 
of  the  side  hulls.  The  troughing  situation  is  especially  relevant  in  the  case  where  the  outer  hulls 
are  relatively  short  compared  to  the  main  hull.  The  situation  is  similar  for  a  cresting  wave 
creating  an  upward  vertical  force  on  the  outer  hulls  due  to  increased  immersion  and  buoyancy. 

The  coincident  side  hull  troughing  and  cresting  conditions  arise  from  the  configuration  of 
the  trimaran’s  side  hulls  not  being  the  same  length  as  the  main  hull.  Depending  on  the  length  of 
the  main  and  side  hulls,  the  side  hulls  could  experience  radical  changes  in  their  buoyant  support 
and  contribution  while  the  main  hull  is  relatively  unaffected  from  the  view  point  of  a  traditional 
primary  bending  conditions. 

3.1.3  Side  Hull  Torsion  in  Longitudinal  Waves 

As  with  the  case  of  side  hull  troughing  and  cresting,  the  unequal  lengths  of  the  hulls 
introduces  a  torsional  structural  loading  on  the  side  hull  and  cross-deck  connection  which  is  not 


experienced  in  traditional  monohull  ship  designs.  The  side  hull  torsional  loading  can  be  seen 
below  in  Figure  6  through  Figure  8. 


Figure  6  -  Negative  Phase  Twisting  with  Side  Hulls  Aft 


Figure  7  -  Positive  Phase  Twisting  with  Side  Hulls  Aft 

Figure  6  through  Figure  8  show  a  profile  view  of  various  possible  configurations  of 
longitudinal  torsional  loading  for  a  trimaran  of  typical  proportions  as  determined  from  Table  2. 
The  dotted  lines  show  the  still  water  position  of  the  trimaran  and  waterline,  the  dashed  line 
shows  the  wave,  and  the  solid  lines  show  the  trimaran’s  response  to  the  wave. 

The  longitudinal  side  hull  torsion  loading  arises  from  similar  conditions  as  that  of  side 
hull  troughing  and  sagging.  For  certain  side  hull  lengths  and  corresponding  wavelengths 
produced  by  a  particular  sea-state,  the  side  hull  and  cross-deck  structure  will  experience  large 
torsional  forces  while  the  main  hull  is  relatively  structurally  unaffected  when  compared  to  its 
worst  case  loading  conditions. 
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Figure  8  -  Twisting  with  Side  Hulls  Amidships 


A  traditional  monohull  ship’s  structural  loading  is  primarily  characterized  by  the 
longitudinal  bending  cases  of  sagging  and  hogging.  However,  the  form  of  the  trimaran  with  its 
outer  stabilizing  hulls  necessarily  involves  two  new  cases  of  structural  loading  in  transverse 
bending  that  is  not  a  concern  with  a  traditional  monohull  design.  These  could  probably  be  best 
described  as  “transverse  hogging”  and  “transverse  sagging.”  The  transverse  sagging  load  case  is 
shown  below  in  Figure  9,  while  transverse  hogging  is  shown  in  Figure  10.  Figure  9  and  Figure 
10  show  an  end  view  of  a  typical  trimaran  ship  represented  with  box  hulls. 


Figure  9  -  Transverse  Sagging 

The  dotted  lines  in  the  figures  show  the  original  still  water  position  of  the  trimaran,  and 
the  heavy  lines  show  the  effect  on  the  ship  in  a  transverse  beam  wave  condition.  Like  the 
loading  conditions  previously  mentioned,  the  transverse  beam  wave  condition  gives  rise  to  a 
situation  where  the  side  hulls  are  again  gaining  or  losing  buoyancy  as  compared  to  the  still  water 
case.  This  change  in  buoyancy  leads  to  vertical  forces  applied  to  the  trimaran’s  cross-deck 
structure.  Depending  on  the  trimaran’s  outer  hull  placement  and  the  wavelengths  of  the  beam 
waves  encountered,  the  side  hulls  could  be  completely  broached  from  or  immersed  in  the  water. 
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Figure  10  -  Transverse  Hogging 

The  “transverse  hogging”  case  in  Figure  10  is  in  reality  an  unstable  condition,  since  the 
trimaran  in  this  loading  condition  will  most  likely  roll  until  one  of  its  outer  hulls  rests  in  the 
water.  However,  the  transverse  hogging  has  been  presented  in  its  unstable  condition  because  it  is 
equivalent  to  the  limiting  load  case  where  the  outer  hull  is  completely  broached  from  the  water. 

3.1.5  Other  Possible  Loading  Cases 

The  structural  loading  conditions  outlined  in  sections  3.1.2,  3.1.3,  and  3.1.4  above  only 
take  into  account  the  limited  scenarios  where  the  wave  front  direction  is  either  completely  in  a 
head  seas  or  beam  seas  orientation.  Obviously,  the  majority  of  the  wave  fronts  experienced  by 
any  ship  in  its  service  life  time  will  in  general  be  at  some  oblique  angle  to  the  ship’s  heading. 
This  oblique  angle  will  contain  some  of  the  characteristics  of  the  load  cases  outline  above  with 
additional  hydrodynamic  forces  tending  to  squeeze  and  spread  the  hulls  in  the  transverse 
direction. 

In  addition  to  the  transverse  hydrodynamic  forces,  other  dynamic  loading  effects  due  to 
the  ship’s  speed  and  heading  will  also  need  to  be  accounted  for  to  perform  a  complete  analysis. 
While  these  other  structural  loadings  and  effects  are  important  to  the  final  structural  design  and 
integrity  of  the  completed  trimaran  ship,  the  complexity  involved  in  accounting  for  each  effect  is 
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Chapter  4  Analysis  Tools 


4.1  Analytical  Trimaran  Model 

Over  the  course  of  this  research  almost  40  versions  of  analytical  trimaran  models  were 
developed  with  each  version  adding  successively  more  functionality  and  detail.  However,  the 
core  of  the  analytical  thesis  model  is  that  it  approximates  a  trimaran  as  three  rigidly  connected 
“box  barges”  of  appropriate  dimensions  obtained  from  Table  2.  The  analytical  model  then 
statically  balances  the  trimaran  on  a  wave  and  calculates  the  load  forces,  as  described  in  section 
3.1,  applied  to  the  cross-deck  structure  between  the  hulls.  In  this  section  the  details  of  the  final 
analytical  thesis  model  will  be  discussed.  Specific  equations  from  the  model  will  be  provided  in 
OAppendix  A.  While  any  programming  language  would  have  been  perfectly  acceptable 
computational  tool,  the  MathCAD  program  was  used  to  perform  the  calculations  for  the 
analytical  model  due  to  its  visually  attractive  mathematical  interface. 

Classification  societies  such  as  the  American  Bureau  of  Shipping  (ABS)  [3]  or  Det 
Norske  Veritas  (DNV)  [4]  often  state  design  codes  for  ships  strength  in  terms  of  the  ship’s 
relevant  dimensions  and  a  sea-state  coefficient  dependent  on  ships  dimensions  to  provide 
sufficient  design  margin  for  the  life  of  a  ship.  The  goal  of  the  analytical  thesis  model  is  to  be 
able  to  state  similar  design  loadings  based  on  the  relevant  parameters  of  a  trimaran  ship.  The 
design  loadings  obtained  for  the  load  cases  described  sections  3.1.2, 3.1.3,  and  3.1.4  are  shown 
below  in  Table  3  while  the  relevant  design  parameters  or  the  trimaran  ship  are  shown  in  Table  4. 


Table  3  -Structural  Loading  from  sections  3.1.2, 3.1.3,  and  3.1.4 


MEM 

the  vertical  force  on  cross-structure  in  longitudinal  waves  (troughing  and  cresting) 

ten 

the  moment  on  cross-structure  in  longitudinal  waves  (longitudinal  twisting) 

re™ 

the  vertical  force  in  transverse  waves  (transverse  hogging/sagging) 

Table  4  -  Relevant  Trimaran  Parameters  Affecting  Design  Loading 


Bm 

the  main  hull’s  beam 

Lm 

the  main  hull’s  length 

Tm 

the  main  hull’s  draft 

Fm 

the  main  hull’s  freeboard  (also  side  hull  freeboard) 

Bs 

the  side  hull’s  beam 

Ls 

the  side  hull’s  length 

the  side  hull’s  draft 

xs 

the  longitudinal  position  of  the  side  hull  wrt  main  hull  amidships 

Ys 

the  transverse  position  of  the  side  hull  wrt  main  hull  centerline 

4.1.1  Analytical  Model  Parts 


The  analytical  model  used  for  this  work  consists  of  two  major  components.  These  are  the 
calculation  of  the  ships  motions  and  forces  and  the  calculation  of  the  representative  sea  states. 
Each  part  is  described  below  with  its  corresponding  equation  provided  in  Appendix  A. 

4. 1.1.1  Entry  of  Statistical  Sea  State 

A  joint  frequency  table  in  both  wave  period  and  significant  wave  height  for  the  northern 
North  Atlantic  Ocean  was  obtained  from  [14],  This  table  was  converted  to  a  joint  probably  table 
whose  results  are  shown  below  in  Table  5  and  Table  6.  This  joint  probability  table  of  wave 
height  and  period  was  used  for  a  statistical  representation  of  the  most  severe  sea  state  in  which 
the  trimaran  hull  form  would  be  subjected.  The  data  represented  in  Table  5  and  Table  6  is  often 
approximated  with  well  known  distributions  such  as  the  Brentschneider  or  other  spectra. 
However,  these  spectra  often  described  in  only  a  few  parameters  and  do  not  always  fully  capture 
the  true  nature  of  the  joint  probability  of  the  sea  state  in  the  wave  height  and  wavelength 
parameters.  Therefore,  the  observed  tabular  data  was  utilized. 

Table  5  -  Joint  Probability  Mass  Table  of  Significant  Wave  Height  and  Period 

Spectral  Peak  Period  (s) 


3 

4 

5 

6 

7 

8 

. 9 

10 

11 

0.5 

0.0006 

0.004 

0.0106 

0.0157 

0.0163 

0.0136 

0.0098 

0.0064 

0.0039 

1.5 

9E-05 

0.0021 

0.0123 

0.0322 

0.0511 

0.0581 

0.0528 

0.041 

0.0285 

2.5 

0 

8E-05 

0.0015 

0.0083 

0.0229 

0.039 

0.0471 

0.0446 

0.0353 

3.5 

0 

0 

6E-05 

0.0008 

0.0048 

0.0137 

0.0241 

0.0296 

0.028 

4.5 

0 

0 

0 

4E-05 

0.0006 

0.0031 

0.009 

0.0156 

0.0188 

5.5 

0 

0 

0 

0 

3E-05 

0.0004 

0.0021 

0.0057 

0.0095 

6.5 

0 

0 

0 

0 

0 

2E-05 

0.0003 

0.0014 

0.0035 

7.5 

0 

0 

0 

0 

0 

0 

2E-05 

0.0002 

0.0009 

8.5 

mm 

0 

0 

0 

0 

0 

0 

2E-05 

EHO 

9.5 

0 

0 

0 

0 

0 

0 

0 

0 

2E-05 

10.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

11.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

12.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

13.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

14.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Where  each  wave  height  indicated  is  the  center  of  a  1  meter  range  of  heights 


Table  6  -  Joint  Probability  Mass  Table  of  Significant  Wave  Height  and  Period 


m 


♦ 


+ 


Spectral  Peak  Period  (s) 


12 

13 

14 

15 

16 

17 

18 

19 

21 

22 

0.5 

0.0023 

0.0013 

0.0007 

0.0004 

0.0002 

0.0001 

7E-05 

4E-05 

2E-05 

2E-05 

1.5 

0.0182 

0.011 

0.0063 

0.0035 

0.0019 

0.001 

0.0006 

0.0003 

0.0002 

0.0002 

2.5 

0.0245 

0.0154 

0.009 

0.005 

0.0026 

0.0013 

0.0007 

0.0003 

0.0002 

0.0001 

3.5 

0.0216 

0.0144 

0.0085 

0.0046 

0.0023 

0.0011 

0.0005 

0.0002 

IE-04 

7E-05 

EBI 

0.017 

0.0123 

0.0075 

0.004 

0.0019 

0.0008 

0.0003 

0.0001 

5E-05 

3E-05 

5.5 

0.0107 

0.0088 

0.0057 

0.0031 

0.0014 

0.0006 

0.0002 

7E-05 

2E-05 

IE-05 

6.5 

0.0053 

0.0053 

0.0039 

0.0022 

0.001 

0.0004 

0.0001 

4E-05 

IE-05 

0 

1151 

0.002 

0.0026 

0.0023 

0.0014 

0.0006 

0.0002 

7E-05 

2E-05 

0 

0 

8.5 

0.0005 

0.001 

0.0011 

0.0008 

0.0004 

0.0001 

4E-05 

IE-05 

0 

0 

9.5 

0.0001 

0.0003 

0.0004 

0.0004 

0.0002 

8E-05 

2E-05 

IE-05 

0 

0 

10.5 

2E-05 

7E-05 

0.0001 

0.0002 

0.0001 

5E-05 

IE-05 

0 

0 

0 

11.5 

0 

IE-05 

4E-05 

6E-05 

5E-05 

2E-05 

IE-05 

0 

0 

0 

12.5 

0 

0 

IE-05 

2E-05 

2E-05 

IE-05 

0 

0 

0 

0 

13.5 

0 

0 

0 

0 

IE-05 

0 

0 

0 

0 

0 

14.5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Where  each  wave  height  indicated  is  the  center  of  a  1  meter  range  of  heights 


4. 1.1. 2  Conversion  of  Statistical  Sea  State 

The  data  represented  in  Table  5  and  Table  6  give  the  joint  probability  of  a  sea  state 
having  a  given  significant  wave  height  and  period.  However,  the  data  given  is  for  short  term 
observations  of  a  narrow-banded,  fully  developed  sea  state  reported  in  terms  of  significant  wave 
height.  However,  significant  wave  height  is  a  one  parameter  descriptor  of  the  probability 
distribution  for  a  fully  developed  sea  state.  The  probability  of  the  peak  values  of  the  actual 
wave  amplitudes  for  a  fully  developed  sea  state  are  described  by  a  Rayleigh  distribution  with 


significant  wave  height  as  the  distribution  function  parameter.  Correspondingly,  each  entry  of 
Table  5  and  Table  6  actually  describes  the  joint  probability  of  a  certain  Rayleigh  distribution 
occurring.  Since  the  limiting  design  of  the  trimaran  is  concerned  with  an  overall  probability  of 


the  wave  loading  conditions  encountered,  the  data  in  Table  5  and  Table  6  was  converted  to 


represent  an  absolute  joint  probability  of  actual  wave  height  verses  wave  period  in  lieu  of 
significant  wave  height. 
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Using  [15],  equation  (1)  was  derived.  Equation  (1)  is  the  Rayleigh  probability 
distribution  of  actual  wave  height  in  terms  of  significant  wave  height. 

P  ray(h  w>1)  =  4' - “T 

(H°.) 

a) 

Where: 

pray  -  the  Rayleigh  distribution  probability  as  functions  of  wave  height  and  the  index  t 

hw  -  the  actual  wave  height 

Hs  -  the  observed  significant  wave  height 

t  -  the  index  of  the  observed  significant  wave  height  in  Table  5  and  Table  6 

Using  the  Rayleigh  probability  distribution  in  equation  (1),  for  each  significant  observed 
wave  height,  the  probability  of  being  in  each  wave  height  range  of  Table  5  and  Table  6  was 
calculated  by  integrating  the  Rayleigh  distribution  in  one  meter  segments.  The  result  of  the 
integral  was  multiplied  by  the  joint  probability  shown  in  Table  5  and  Table  6  and  then  added  to 
the  other  associated  probabilities  affecting  that  wave  height  region.  The  result  is  shown  in  below 
in  Table  7  and  Table  8. 


Table  7  -  Joint  Probability  of  True  Wave  Height  and  Period  (3-11  second  periods) 

Spectral  Peak  Period  (s) 


5 

6 

0.0183 

0.0371 

0.0054 

0.0163 

0.0007 

0.003 

9E-05 

0.0006 

IE-05 

IE-04 

IE-06 

2E-05 

IE-07 

3E-06 

2E-08 

5E-07 

0 

6E-08 

0 

IE-08 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

IE-08 


4E-09 


0 


3E-08  5E-07  4E-06 


8E-09  IE-07  IE-06 


2E-08  3E-07 


5E-09  IE-07 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


10 


0.0492 


0.0502 


0.0252 


0.0115 


0.005 


0.0021 


0.0008 


0.0003 


0.0001 


5E-05 


2E-05 


6E-06 


2E-06 


7E-07 


2E-07 


7E-08 


IE-08 


4E-09 


11 


0.0372 


0.0424 


0.0245 


0.0128 


0.0063 


0.003 


0.0014 


0.0006 


0.0003 


0.0001 


5E-05 


2E-05 


7E-06 


3E-06 


IE-06 


4E-07 


IE-07 


5E-08 


9E-09 


4E-09 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 
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Table  8  -  Joint  Probability  of  True  Wave  Height  and  Period  (12-22  second  periods' 

Spectra!  Peak  Period  (s 


13 

14 

15 

16 

17 

18 

19 

0.0163 

0.0096 

0.0053 

0.0028 

0,0014 

0.0007 

0.0004 

0.0212 

0.0128 

0.007 

0.0036 

0.0018 

0.0008 

0.0004 

0.0145 

0.009 

0.005 

0.0025 

0.001 1 

0.0005 

0.0002 

0.0089 

0.0058 

0.0032 

0.0016 

0.0007 

0.0003 

0.0001 

0.0052 

0.0035 

0.002 

0.001 

0.0004 

0.0001 

5E-05 

0.003 

0.0021 

0.0012 

0.0006 

0.0002 

8E-05 

3E-05 

0.0016 

0.0012 

0.0007 

0.0004 

0.0001 

5E-05 

IE-05 

0.0009 

0.0007 

0.0004 

0.0002 

8E-05 

3E-05 

8E-06 

0.0005 

0.0004 

0.0003 

0.0001 

5E-05 

IE-05 

4E-06 

0.0002 

0.0002 

0.0002 

8E-05 

3E-05 

9E-06 

2E-06 

0.0001 

0.0001 

9E-05 

5E-05 

2E-05 

5E-06 

IE-06 

6E-05 

7E-05 

5E-05 

3E-05 

IE-05 

3E-06 

7E-07 

3E-05 

3E-05 

3E-05 

2E-05 

6E-06 

2E-06 

4E-07 

IE-05 

2E-05 

2E-05 

IE-05 

4E-06 

IE-06 

2E-07 

7E-06 

IE-05 

9E-06 

6E-06 

2E-06 

hii:g»M 

IE-07 

3E-06 

5E-06 

5E-06 

3E-06 

IE-06 

3E-07 

5E-08 

IE-06 

2E-06 

3E-06 

2E-06 

7E-07 

2E-07 

2E-08 

7E-07 

IE-06 

IE-06 

IE-06 

4E-07 

IE-07 

IE-08 

3E-07 

6E-07 

7E-07 

6E-07 

2E-07 

5E-08 

4E-09 

IE-07 

3E-07 

4E-07 

4E-07 

IE-07 

3E-08 

2E-09 

4E-08 

IE-07 

2E-07 

2E-07 

6E-08 

IE-08 

0 

2E-08 

7E-08 

IE-07 

IE-07 

4E-08 

8E-09 

0 

3E-09 

2E-08 

4E-08 

5E-08 

2E-08 

3E-09 

0 

2E-09 

IE-08 

2E-08 

3E-08 

9E-09 

2E-09 

0 

0 

3E-09 

6E-09 

IE-08 

3E-09 

0 

0 

0 

2E-09 

3E-09 

8E-09 

2E-09 

0 

0 

0 

0 

0 

3E-09 

0 

0 

0 

0 

0 

0 

2E-09 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

IE-06 


5E-07 


2E-07 


8E-08 


3E-08 


8E-09 


2E-09 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


IE-07 


IE-08 


3E-09 


0 


0 


0 


0 


4.1. 1.3  Equations  of  Applied  Waves 

To  create  the  loading  conditions  discussed  above  in  sections  3.1.2, 3,1.3,  and  3.1.4  from 
the  wave  height  and  period  data  in  Table  7  and  Table  8,  the  waves  were  applied  as  described 
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First  the  deep-water  wavelength  approximation  was  made  by  equation  (2)  for  each  period 
in  Table  7  and  Table  8. 

A=JL.(t)2 

2-71 

(2) 

Where: 

X  -  the  wavelength  of  the  wave 

g  -  the  gravitational  constant 

T  -  the  wave  period  from  Table  7  and  Table  8 

Equation  (2)  applied  to  each  indicated  wave  period  yielded  the  wavelengths  indicated  in 
below  Table  9. 


Table  9  -  Deep  Water  Wavelengths 


Period 

(s) 

Wavelength 

(m) 

3 

14.05 

4 

24.97 

5 

39.02 

6 

56.19 

7 

76.48 

8 

99.89 

9 

126.42 

10 

156.08 

11 

188.85 

12 

224.75 

13 

263.77 

14 

305.91 

15 

351.17 

16 

399.56 

17 

451.06 

18 

505.69 

19 

563.44 

21 

688.30 

22 

755.42 

Next  the  following  wave  equations  for  the  longitudinal  and  transverse  directions  were 
applied  with  the  proper  phasing  to  apply  the  loading  conditions  described  in  sections  3.1.2,  3.1.3, 
and  3.1.4.  The  coordinate  system  used  for  the  analytical  model  was  those  of  traditional  ship 
coordinate  systems  of  x-longitudinal,  y-transverse,  and  z-vertical.  The  coordinate  system  of  the 
analytical  model  is  shown  below  in  Figure  11. 


z  -  axis 


Figure  11  -  Analytical  Model  Coordinate  System 


The  wave  equations  used  to  create  the  specific  loading  conditions  are  show  below  in 
equations  (3)  and  (4). 


f  °  1 

'"posjwist  N 

^  wave_lg(x’P’<J»<l,»  -^s 

^  Wn  f0.1T 

)  = - ^  sin - (x-  X„)  -  <}> 

1  2  |_Xq  V  S' 

II 

-&• 

90 

180 

,270, 

•deg 

^lg11 

trough 

neg_twist 

,  crest  j 

(3) 


w_ 


^  wave_tr^y’P ’4  j^)  •  ^  cos 


~(y) 

LAq 


0  ]  ( hog 

♦»  =  l  180  rg  ♦‘'Hsag 


Where: 

Hwavejg  -  is  the  longitudinal  wave  amplitude  as  functions  of  the  shown  variables 
HwaVe_tr  -  is  the  transverse  wave  amplitude  as  functions  of  the  shown  variables 
hw  -  the  actual  wave  height  from  Table  7  and  Table  8 
x  -  the  longitudinal  position 

Xs  -  the  longitudinal  position  of  the  side  hull  wrt  main  hull  amidships 
y  -  the  transverse  position 
k  -  the  wavelength  of  the  wave 

cp  -  the  phase  of  the  wave  applied  for  each  loading  condition 
p  —  the  index  of  the  wave  height  Table  7  and  Table  8 
q  -  the  index  of  the  wavelength  from  Table  9 


Equations  (3)  and  (4)  define  the  wave  amplitude  in  the  longitudinal  and  transverse 
directions  respectively.  Both  equations  are  functions  of  the  index  p  which  references  the  wave 


heights  in  Table  7  and  Table  8.  The  wave  definition  equations  are  also  functions  of  the  index  q 
referencing  the  wavelengths  in  Table  9.  Additionally  equation  (3)  is  defined  so  that  the  origin  of 
the  wave  is  longitudinally  offset  to  coincide  with  the  longitudinal  center  of  the  side  hull,  Xs. 
Equations  (3)  and  (4)  are  also  defined  as  a  function  of  applied  phase  angle  so  that  each  phase 
angle  can  be  easily  applied  to  generate  the  load  cases  defined  in  sections  3.1.2,  3.1.3,  and  3.1.4. 
Finally,  sinusoidal  waves  were  used  to  approximate  the  wave  functions  rather  than  trochoidal 
waves,  which  more  truly  represent  water  waves,  since  the  difference  between  sinusoidal  and 
trochoidal  waves  for  the  wavelengths  of  concern  is  relatively  small. 


4. 1.1. 4  Calculation  of  Motion  due  to  Applied  Waves 

After  the  waves  applied  to  the  trimaran  hull  form  were  defined  as  in  equations  (3)  and  (4) 
above,  the  next  step  in  finding  the  load  forces  on  the  trimarans  cross-structure  was  to  perform  a 
static  balance  of  the  ship  for  a  particular  wave.  For  the  longitudinal  wave  cases,  the  static 
balance  required  simultaneously  solving  for  the  pitch  and  heave  of  the  ship.  While  solving 
simultaneously  equations  for  heave  and  pitch  of  a  ship  is  conceptually  simple,  the  practice  of 
executing  this  concept  is  somewhat  cumbersome.  To  perform  the  static  balance,  the  underwater 
volume  and  the  first  moment  of  that  volume  subjected  to  the  wave  must  balance  with  the 
displacement  and  moment  of  the  ship  in  the  pitched  and  heaved  condition.  This  process  involved 
iteratively  solving  equations  (5)  and  (6)  below  using  a  modified  Newton’s  method.  Equation  (5) 
is  the  heave  equation. 


v  still  -  B  m 


wavejg^^  AT  (x  X  cf)sin(©  pitch)]  dx  — 


wavejgta*^  AT  (x  X  cf)-sin|©  pjtch)J  dx 


Where: 

AT  -  the  change  in  heave  of  the  ship  due  to  the  longitudinal  wave 

©pitch  -  the  change  in  pitch  of  the  ship  due  to  the  longitudinal  wave 

HWavejg  -  is  the  longitudinal  wave  height  as  functions  of  position  and  phase  from  (3)* 

x  -  the  longitudinal  position 

tp  -  the  phase  of  the  wave  applied  for  each  loading  condition 
Bm  -  the  main  hull’s  beam 
Lm  -  the  main  hull’s  length 
Bs  -  the  side  hull’s  beam 
Ls  -  the  side  hull’s  length 

Xs  -  the  longitudinal  position  of  the  side  hull  wrt  main  hull  amidships 
Xcf-  the  longitudinal  position  of  the  entire  ship’s  center  of  floatation 
Vsti„  -  the  original  still  water  displacement 
Trinistiii  -  the  original  still  water  trim  (assumed  to  be  zero) 

*  For  simplicity  the  wave  height  here  is  only  shown  as  functions  of  position  and  phase  with  wave  height 
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Equation  (6)  is  the  heave  equation. 


Trimstill=  B  m 


2 

-L. 


[H  wavejg^  "  AT  -  (x  -  X  cf)sin(©  pitch)]  (x  -  X  cf)  dx . 


+  2B, 


rxs+- 


X, 


[H  wave_lg(x><l>)  "  AT  -  (x  -  Xcf).sin(©  pitch)]  (x  -  Xcf)  dx 


Where: 

The  definitions  of  terms  from  equation  (5)  apply. 


(6) 


As  with  the  longitudinal  case  above,  the  heave  of  the  trimaran  due  to  transverse  waves 
was  calculated  in  a  similar  manner.  However,  since  this  work  does  not  deal  with  the  affects  of 
roll  on  the  ship  and  takes  as  a  worst  case  the  “transverse  hogging  and  sagging”  described  in 
section  3.1.4,  only  one  integral  equation  was  needed  to  solve  for  the  ship’s  heave  in  transverse 
waves. 

,.B_i 

^  still  =  Lm' 


(Hwave_tr(y^)-AT)dy  +  2-Ls- 


-B, 


(HwaveJrM  ~AT)dy 


•'y 


(7) 

Where: 

AT  -  the  change  in  heave  of  the  ship  due  to  the  transverse  wave 

Hwave_tt-  is  the  transverse  wave  height  as  functions  of  position  and  phase  from  (4)* 

y  -  the  transverse  position 

<p  -  the  phase  of  the  wave  applied  for  each  loading  condition 
Bm  -  the  main  hull’s  beam 
Lm  -  the  main  hull’s  length 
Bs  -  the  side  hull’s  beam 
Ls  -  the  side  hull’s  length 

Ys  -  the  transverse  position  of  the  side  hull  wrt  main  hull  centerline 
VstiU  -  the  original  still  water  displacement 

*  For  simplicity  the  wave  height  here  is  only  shown  as  functions  of  position  and  phase  with  wave  height 
and  wavelength  indices  omitted. 


While  equations  (5),  (6),  and  (7)  display  the  conditions  for  a  static  balance  on  waves,  the 
motion  parameters  of  heave  and  pitch  (AT  and  ©pitch)  were  solved  for  in  terms  of  a  change  from 
the  initial  still  water  conditions  that  were  assumed  to  be  zero.  The  motion  parameters  of  heave 
and  pitch  in  the  analytical  model  were  left  as  functions  of  several  variables  that  defined  the 
general  characteristics  of  the  trimarans  dimensions,  so  that  the  motion  of  several  different 
variations  of  trimarans  could  be  computed. 


4. 1 . 1 .5  Calculation  of  Forces  due  to  Motion  in  W aves 

Once  the  motions  of  the  trimaran  due  to  the  static  balance  were  found  as  explained 
section  4. 1.1. 4,  the  forces  due  to  those  motions  were  calculated.  The  motion  parameters  of  the 
ship  in  waves  in  section  4. 1.1. 4  were  determined  in  reference  to  the  still  water  condition. 
Consequently,  the  force  calculation  was  performed  as  a  change  from  the  still  water  waterline 
position  of  the  trimaran’s  side  hulls.  For  example,  if  the  side  hull  is  more  immersed  than  the  still 
water  position  as  in  the  “transverse  sagging”  case  shown  in  Figure  12,  then  the  buoyancy  of  the 
water  provides  an  upward  force  on  the  cross-structure. 


Figure  12  -  Transverse  Sagging  Applying  Upward  Buoyant  Force  on  Cross  Structure 
The  force  calculations  for  transverse  hogging  and  longitudinal  cresting  or  troughing  are 
similar.  In  these  cases  the  particular  wave  height  and  phase  that  the  trimaran  encounters  either 
provides  more  or  less  buoyancy  to  the  side  hulls  and  consequently  cantilevered  type  vertical 
forces  are  applied  to  the  cross  structure. 

The  twisting  force  applied  to  the  cross  structure  of  the  trimaran  was  calculated  in  a 
similar  manner  as  the  vertical  forces.  The  ship  was  statically  balanced  on  the  wave  at  phases  of 
maximum  twisting  forces  on  the  cross  structure,  and  then  the  twisting  moment  due  to  the  wave 
was  calculated  as  a  change  from  the  still  water  condition.  A  case  of  this  twisting  situation  is 
shown  in  Figure  13  which  causes  a  positive  twisting  moment  to  be  applied  to  the  cross  structure. 


The  position  of  the  calculation  of  the  twisting  moment  coincided  with  the  Ion 
center  of  the  side  hull.  The  vertical  position  of  the  twisting  moment  was  assumed  to 
original  vertical  position  of  side  hull  draft,  which  is  an  approximate  position  due  the 


Figure  13  -  Positive  Phase  Twisting  with  Side  Hulls  Aft 

The  analytical  model  calculated  the  vertical  forces  and  twisting  forces  applied  to  the  side 
hull  due  to  the  motion  in  waves.  The  model  was  generalized  such  that  these  forces  were 
calculated  as  functions  of  many  variables  so  that  the  model  could  be  used  to  repetitively 
calculated  forces  for  different  wave  states  and  trimaran  configurations.  These  parameters 
included  those  of  wave  height,  wavelength,  and  wave  phase  from  the  wave  equations  of  (3)  and 
(4).  The  pertinent  parameters  of  the  trimaran  such  as  length,  beam,  draft,  hull  spacing  and 
freeboard  were  also  included  in  the  force  calculations. 
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The  dependent  forces  variables  and  there  independent  variables  are  shown  below  in 


equations  (8),  (9),  and  (10)  with  position  of  application  shown  in  Figure  14  their  full  definition  of 


calculation  shown  in  Appendix  A. 

FBlg(p)qAXs,Xcf,Ts,Ls,Bs,Fm,Tm,Bm,Ln]) 


(8) 


MB  lg(p>4>$>  Xs,Xcf,T  S,LS,B 


(9) 


Where: 

FBig  -  the  vertical  force  on  cross-structure  in  longitudinal  waves  (troughing  and  cresting) 
MBlg  -  the  moment  on  cross-structure  in  longitudinal  waves  (longitudinal  twisting) 

FBtr  -  the  vertical  force  in  transverse  waves  (transverse  hogging/sagging) 

p  -  the  index  of  the  wave  height  Table  7  and  Table  8 

q  -  the  index  of  the  wavelength  from  Table  9 

cp  -  the  phase  of  the  wave  applied  for  each  loading  condition 

X,.f  -  the  longitudinal  position  of  the  entire  ship’s  center  of  floatation 

Xs  -  the  longitudinal  position  of  the  side  hull  wrt  main  hull  amidships 

Tm  -  the  main  hull’s  draft 

Bm  -  the  main  hull’s  beam 

Lm  -  the  main  hull’s  length 

Fm  -  the  main  hull’s  freeboard  (also  side  hull  freeboard) 

Ts  -  the  side  hull’s  draft 
Bs  -  the  side  hull’s  beam 
Ls  -  the  side  hull’s  length 

Ys  -  the  transverse  position  of  the  side  hull  wrt  main  hull  centerline 


Figure  14  -  Analytical  Model  Position  of  Calculated  Forces 


One  main  assumption  for  the  calculation  of  forces  and  moments  in  the  analytical  model 
was  that  the  initial  still  water  position  of  the  side  hulls  provided  a  neutrally  buoyant 
displacement.  This  neutrally  buoyant  condition  enforces  a  condition  that  there  are  no  applied 
forces  to  the  cross  structure  due  to  still  water  buoyancy  conditions.  While  this  would  seem  to  be 
a  reasonable  design  condition,  it  would  not  be  a  necessary  condition  and  any  imbalance  in  the 
still  water  forces  would  require  accounted  for.  For  the  purpose  of  this  work,  however,  the  still 
water  loading  was  assumed  to  be  zero  so  that  the  nature  of  the  structural  loading  in  waves  could 
be  determined  by  equations  (8),  (9),  and  (10). 

Another  simplification  in  the  analytical  model  was  that  wet  deck  slamming  of  waves  with 
large  wave  heights  into  the  trimaran  cross-deck  structure  was  not  taken  into  account.  While  this 
assumption  will  eventually  need  to  be  restored  to  fully  calculate  the  structural  effects  of  sea  state 
on  trimarans,  the  effects  were  not  found  to  be  extremely  relevant  given  the  scope  of  this 
investigation.  For  example,  for  a  100m  long  trimaran  of  representative  dimensions  from  Table  1, 
to  obtain  impact  of  waves  with  the  trimarans  cross  deck  structure  would  require  significant  wave 
heights  on  the  order  of  6  meters.  This  wave  height  corresponds  to  sea  states  of  6  or  7  [14].  Not 
only  are  these  sea  states  very  high  for  a  ship  of  100  meters,  but  also  these  sea  states  are  relatively 
unlikely  probabilistically  as  can  be  seen  in  Table  5  and  Table  6.  Therefore,  discounting  the 
cross-deck  slamming  for  this  analysis  can  be  justified  for  this  work. 

{  ■  •  ,  if,  '  •  '  .  : 

4.1. 1.6  Calculation  of  Design  Forces  due  to  a  Full  Sea  Spectrum 

Once  the  functions  used  to  calculate  the  forces  on  the  cross  deck  structure  were  defined 
in  equations  (8),  (9),  and  (10),  it  was  then  possible  to  define  a  design  load  that  the  cross  deck 
structure  of  a  trimaran  must  withstand.  Design  loads  are  often  stated  in  terms  of  safety  factors  or 
design  margins.  But  at  the  root  of  a  statistical  design  process  design  loads  can  also  be  stated  in 
terms  of  the  reliability  index,  which  is  what  was  used  in  this  work.  The  reliability  index  used  for 
this  investigation  was  5,  which  is  suggested  for  naval  war  ships  in  [16].  In  the  strictest  sense,  the 
reliability  index  must  take  into  account  the  probability  distributions  of  both  the  load  and  the 
strength  factors  of  a  design.  However,  for  the  purpose  of  this  work  only  the  probability 
distribution  of  the  sea-state  and  thus  wave  loading  was  used  in  the  calculation  of  the  reliability 
index.  Omitting  the  probabilistic  nature  of  the  strength  curves  equates  to  knowing  the  failure 
strength  deterministically,  which  is  a  reasonable  simplification  when  considering  steel 


manufacturing  quality  assurance  techniques.  Once  the  design  loading  for  a  trimaran  of  general 
dimensions  is  determined,  the  partial  safety  factors  for  each  possible  failure  mode  could  then  be 
calculated  if  desired  [16]. 

In  addition  to  the  reliability  index,  the  other  information  need  to  determine  the  design 
loading  for  each  load  case  of  defined  in  sections  3.1.2, 3.1.3,  and  3.1.4,  was  the  susceptible 
wavelength  for  trimaran  of  particular  dimensions.  For  example,  the  United  States  Naval  design 
standard  [17]  and  [18]  for  primary  hull  bending  of  a  monohull  ship  has  been  to  design  for  a 
trochoidal  wave  as  shown  in  (11). 


hw=W^ 


X  =  L 


m 


(ID 


While  equation  (11)  approximates  the  loadings  of  a  fully  statistical  sea-state  relatively 
well  for  traditional  naval  ships,  the  formulation  does  not  work  well  for  a  trimaran  due  to  the 
varying  lengths  and  spacing  of  the  trimaran  hulls.  Therefore,  the  susceptible  wave  length  for 
each  load  case  described  in  sections  3.1.2,  3.1.3,  and  3.1.4  was  determined  by  finding  the 
wavelength  in  which  the  maximum  forces  were  produced  from  equations  (8),  (9),  and  (10)  when 
considering  the  fully  statistical  sea  state  as  described  in  Table  7  and  Table  8. 

However,  computing  the  forces  produced  for  each  entry  in  Table  7  and  Table  8  requires 
statically  rebalancing  the  trimaran  for  each  wavelength  and  wave  height  in  Table  7  and  Table  8. 
The  amount  of  computational  effort  to  solve  the  static  motion  balance  for  all  the  entries  in  the 
statistical  sea  state  is  extremely  large  and  essentially  unnecessary  since  the  low  wave  height 
waves  will  not  produce  high  forces  regardless  of  wave  length.  Therefore,  to  determine  the 
wavelengths  to  which  a  trimaran  of  particular  dimensions  were  susceptible,  the  forces 
encountered  for  each  entry  of  a  statistical  sea  state  described  in  Table  7  and  Table  8  were 
effectively  searched  for  the  highest  forces.  To  alleviate  the  need  for  calculating  every  force 
produced  at  every  possible  sea  state  described  in  Table  7  and  Table  8,  two  simplifications  were 
made.  First,  only  wavelengths  that  would  have  the  ability  to  produce  a  maximum  force  were 
considered.  For  example,  the  structural  response  of  a  trimaran  with  a  main  hull  length  of  100 
meters  is  minimal  on  a  wave  of  with  a  wavelength  of  755  meters  from  Table  9.  This  situation 
can  be  seen  in  Figure  15.  In  the  case  where  the  wavelength  encountered  is  much  longer  than  the 
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length  of  the  ship  the  ship  gently  rides  the  wave  without  experiencing  large  structural  loadings. 
This  concept  is  also  the  basis  for  the  Naval  design  rules  in  (1 1). 


Figure  15  -  Profile  View  of  100m  Trimaran  Balanced  on  a  Wave  of  Long  Wavelength 

The  second  method  utilized  to  reduce  the  computational  effort  of  finding  the  susceptible 
wavelengths  of  a  trimaran  was  to  only  calculate  the  forces  produced  from  the  maximum  wave 
height  for  each  respective  wavelength  which  had  a  non-zero  statistical  probability  in  Table  7  and 
Table  8.  Once  the  forces  for  the  maximum  wave  height  for  each  wavelength  were  calculated,  the 
forces  were  compared  and  the  susceptible  wavelength  was  determined  from  the  corresponding 
maximum  force.  For  example,  for  the  wave  periods  of  4  and  5  seconds  the  forces  produced  for 
5.5  and  7.5  meter  wave  heights  were  calculated  and  compared. 

Using  these  methods,  the  forces  for  each  loading  case  in  sections  3.1.2, 3.1.3,  and  3.1.4 
were  calculated  when  subjected  to  applicable  portions  of  the  statistical  sea  state  in  Table  7  and 
Table  8.  From  those  calculations,  the  susceptible  wavelength  index  was  set  to  the  wavelength 
for  the  corresponding  maximum  force  obtained.  Once  the  susceptible  wavelength  for  a  given 
loading  condition  was  determined,  the  design  loading  for  that  mode  of  loading  could  be 
calculated. 

To  determine  the  design  loading  for  the  susceptible  wavelength,  the  joint  probabilities  in 
Table  7  and  Table  8  were  converted  to  a  marginal  probability  distribution  at  the  susceptible 
wavelength.  Next,  the  force  developed  at  each  wave  height  in  Table  7  and  Table  8  at  the 
susceptible  wavelength  was  found.  Finally,  using  the  marginal  probability  and  the  developed 
forces  at  one  meter  wave  height  increments,  the  mean  and  standard  deviation  of  the  force 
distribution  at  the  susceptible  wavelength  was  determined.  From  there  the  design  loading  for 
each  load  case  was  determined  using  the  reliability  index  of  5  so  that  the  design  loading  would 
be  5  standard  deviations  from  the  mean  loading. 

This  methodology  of  determining  design  loading  was  based  on  the  joint  probability  mass 
function  of  Table  7  and  Table  8.  While  this  method  certainly  lends  insight  into  the  structural 
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loading  of  trimarans  due  to  wave  motions,  the  tabular  nature  of  the  probability  distributions  in 
Table  7  and  Table  8  leaves  the  evaluation  of  the  forces  experienced  in  waves  somewhat  granular. 
A  more  in  depth  investigation  of  this  subject  could  be  performed  using  a  continuous  functions  to 
describe  the  sea-state  probability  distribution  as  in  [13],  but  for  the  purposes  of  this  work’s 
investigation  into  the  basic  structural  loading  attributes  of  a  trimaran  hull  for  would  be 
unnecessary. 

4. 1 . 1 .7  Calculation  of  Design  Forces  for  Varying  Sizes  of  Trimarans 

Once  the  design  forces  as  functions  of  the  relevant  parameters  of  the  trimaran  could  be 
determined  using  the  methods  in  section  4. 1.1. 6,  the  final  component  of  the  analytical  thesis 
model  computed  the  design  forces  for  468  different  variations  of  trimaran  sizes  and  placement 
configurations  to  investigate  the  possible  design  space  of  trimaran  hull  forms.  The  inputs  that 
were  allowed  to  vary  were  main  hull  length,  side  hull  length,  side  hull  x-location,  side  hull  y- 
location  these  are  shown  in  Table  10. 


Table  10  -  Parameters  varied  for  Trimaran  Design  Space 


Lm  (main  hull’s  length) 

60  -  300  meters  in  20meter  increments 

Ls  (side  hull’s  length) 

0.25  -  0.50  in  increments  of  0.05  of  L„,  (main 
hull’s  length) 

Xs  (longitudinal  position  of  the  side  hull  with 
respect  to  main  hull  amidships) 

Five  equally  spaced  position  from  aligned 
amidships  to  where  aft  perpendiculars  of  side 
and  main  hulls  are  aligned 

Ys  (transverse  position  of  the  side  hull  with 
respect  to  main  hull  centerline) 

1  -  1.5  of  Bm  (main  hull’s  beam)  in  increments 
of  0.1 

The  other  relevant  parameters  needed  for  the  computations,  for  example  side  hull  beam, 
were  calculated  from  fixed  relations  from  Table  2  and  reference  [1]  these  are  shown  below  in 
Table  11. 


Table  11  -  Parameters  Held  in  Fixed  Relations  for  Trimaran  Design  Space 


Bm  (main  hull’s  beam) 

Maintain  Lm  /  Bm  =  14 

Tm  (main  hull’s  draft) 

Maintain  Bm  /  Tm  =  2 

Dm  (main  hull’s  depth) 

Required  to  find  (Fm ,  main  hull’s  freeboard) 

Maintain  Dm  /  Tm  =  2.4 

Vs  (side  hull  volume) 

Maintain  at  4%  Vm  (main  hull  volume) 

Bs  (side  hull’s  beam) 

Maintain  Bs  =  Ts  (side  hull’s  draft) 

Ts  (side  hull’s  draft) 

Fixed  by  above  constraints  on  Vs  Ls  Bs 

The  results  of  the  required  design  force  calculations  for  the  design  space  are  shown  in 
OAppendix  B  with  an  explanation  of  the  analysis  of  the  data  provided  in  section  5.1. 


4.2  Finite  Element/  MAESTRO  Analysis 

Once  the  results  from  the  analytical  thesis  model  were  obtained,  a  more  refined  model  of 
a  trimaran  was  made.  The  program  used  to  create  a  more  realistic  model  of  a  trimaran  was  the 
ship  structural  design  program  called  Method  for  Analysis  Evaluation  and  Structural 
Optimization  (MAESTRO).  The  MAESTRO  program,  distributed  by  Proteus  Engineering,  was 
developed  to  analyze  ship  structures  in  a  quasi-finite  element  analysis,  with  large  stiffened  panels 
as  the  elements  between  the  finite  element  nodes.  The  utility  of  using  a  ship  structural  program 
such  as  MAESTRO  is  that  it  allows  the  ability  to  automatically  balance  a  floating  ship’s 
structure  in  water  utilizing  linear  wave  theory.  Naturally,  the  ability  to  include  all  the  relevant 
hydrodynamic  forces  in  the  structural  design  of  a  ship  is  critical  to  accurate  and  safe  ship  design. 
In  this  case,  MAESTRO  only  performs  a  static  balance  of  a  ship  in  waves,  rather  than  a  full 
dynamic  analysis  of  the  structure  interacting  with  the  ocean.  However,  allowing  a  static  analysis 
that  includes  corrections  for  linear  wave  theory  makes  MAESTRO’s  results  that  much  closer  to 
reality. 

4.2. 1  Refined  Finite  Element  MAESTRO  Model 

The  complexity  of  building  the  MAESTRO  model  to  the  individual  scantling  of  each 
component  of  the  trimaran  would  be  prohibitive  for  an  initial  investigation  of  the  previous  design 
space  performed  in  4. 1.1. 7.  Therefore,  one  MAESTRO  model  was  made  to  test  the  validity  of 
the  prediction  of  the  analytical  model  presented  in  section  4.1.  The  particulars  of  the  model  in 
MAESTRO  were  obtained  from  [19]  and  are  shown  below  in  Table  12. 

Table  12  -  Particular  Dimensions  of  Real  Ship  Model  Used 


parameter 

meters 

Lm 

106 

Bm 

9 

Tm 

5.2 

Dm 

11.9 

Ls 

35 

Bs 

1.8 

Ts 

1 

Xs 

-35.5  from 
amidships 

Ys 

9 

Using  the  relations  in  the  analytical  thesis  model  developed  in  section  4,1,  the  design 
loadings  for  longitudinal  cresting/troughing,  longitudinal  twisting,  and  transverse  ® 

hogging/sagging  from  sections  3.1.2, 3.1.3,  and  3.1.4,  were  found.  The  results  of  these  design 
values  are  shown  below  in  Table  13. 

Table  13  -  Analytical  Model  Predictions  for  Comparison  to  MAESTRO  0 


Load  Case 

Force  /  Moment 

Susceptible 
Wavelength  (m) 

Longitudinal  Troughing 

-5.49  x  10s  N 

76.5 

Longitudinal  Cresting 

5.72  x  10s  N 

76.5 

Positive  Longitudinal  Twist 

9.12  x  10b  Nm 

56.2 

Negative  Longitudinal  Twist 

-1.16  x  107Nm 

99.9 

Transverse  Sagging 

1.09  x  10°  N 

25.0 

Transverse  Hogging 

-5.22  x  10*  N 

14.0 

A  representative  MAESTRO  model  from  [19]  using  the  parameters  of  Table  12  was 
simulated  using  the  MAESTRO  program.  This  model  is  shown  below  in  Figure  16  through 
Figure  18.  The  input  file  of  the  MAESTRO  model  is  approximately  60  pages  of  text  files,  and 


4.2.2  Statistically  Significant  Waves  for  the  MAESTRO  Model 


As  mentioned  previously  in  section  4. 1.1. 6  a  reliability  index  of  5  was  used  in  the 
calculation  of  the  design  forces  that  a  trimaran’s  cross  structure  should  be  able  to  withstand 
under  the  various  loading  cases  defined  in  sections  3.1.2, 3.1.3,  and  3.1.4.  The  use  of  the 
reliability  index  implies  that  the  ship’s  Response  Amplitude  Operator  (RAO)  to  a  given  sea  state 
is  known  so  that  the  forces  produced  by  a  particular  sea  state  on  the  trimaran  structure  in 

particular  failure  modes  can  be  calculated.  Alternatively,  as  in  the  case  of  the  analytical  model  in 

' 

section  4.1,  the  forces  produced  could  be  calculated  directly  to  response  of  the  fully  statistical 
sea  state  input  without  specifically  knowing  the  RAO.  However,  to  adequately  compare  the 
design  forces  predicted  by  the  analytical  model  to  the  finite  element  model  in  MAESTRO  the 
RAO  of  the  MAESTRO  model  would  have  to  be  known  or  approximated.  However,  calculating 
the  RAO  of  a  new  ship  such  as  the  trimaran  would  require  a  fully  dynamic  analysis  including 
hydrodynamic  effects,  which  is  the  difficult  and  largely  unknown  task  in  the  first  place. 

This  unknown  RAO  for  the  trimaran  leads  to  the  question  of  which  wave  height 
corresponds  to  producing  a  force  that  has  a  reliability  index  of  5,  giving  a  design  force  of  5 
standard  deviations  over  the  mean  force  for  the  fully  statistical  sea  states  described  in  Table  7 
and  Table  8.  When  dealing  with  ship  response  to  waves,  often  motion  RAO’s  are  unity  in  the 
low  frequency  limit  [15],  which  is  exactly  the  static  balance  performed  by  MAESTRO. 

Therefore,  the  assumption  was  made  that  using  wave  heights  that  were  five  standard  deviations 
above  the  mean  wave  height  for  each  wave  length  would  produce  a  force  response  in  the  various 
load  cases  described  in  sections  3.1.2, 3.1.3,  and  3.1.4  that  are  five  standard  deviations  over  the 


Using  the  data  from  Table  7  and  Table  8,  the  wave  height  to  produce  a  response 
equivalent  to  the  reliability  index  of  5  was  calculated  using  the  algorithm  shown  in  Appendix  C 
with  the  results  shown  below  in  Table  14. 


Table  14  -  Equivalent  Wave  Height  for  Reliability  Index  of  5 


Wave  Period 

Wave  Length 

Mean  +  5SD  Wave 

(sec) 

(m) 

Height  (m) 

3 

14.05 

1.81 

4 

24.97 

2.71 

5 

39.02 

3.48 

6 

56.19 

4.24 

7 

76.48 

5.07 

8 

99.89 

5.98 

9 

126.42 

7.00 

10 

156.08 

8.10 

11 

188.85 

9.26 

12 

224.75 

10.40 

13 

263.77 

11.45 

14 

305.91 

12.32 

15 

351.17 

12.85 

16 

399.56 

12.97 

17 

451.06 

12.30 

18 

505.69 

11.10 

19 

563.44 

9.91 

21 

688.30 

8.13 

22 

755.42 

7.08 

The  values  in  Table  14  are  also  plotted  in  Figure  19.  The  shape  and  magnitude  in  Figure 
19  roughly  correspond  to  the  wave  coefficient  values  in  the  rule  based  formulas  for  longitudinal 
bending  in  references  [3]  and  [4],  The  fact  that  Figure  19  is  so  similar  to  current  rule  based 
classification  societies  wave  coefficients  implies  that  the  simplification  of  applying  a  wave  with 
a  wave  height  that  is  5  standard  deviations  over  the  mean  for  a  given  wavelength  is  most  likely  a 
valid  simplification  and  concurs  with  current  classification  society  practices. 


Mean  +  5SD  Wave  Height  (m) 


Wavelength  (m) 

Figure  19  -  Reliability  Index  Wave  Heights  of  5  verses  Wavelength 

Using  the  wave  height  and  length  data  from  Table  14  in  conjunction  with  the  load  cases 
defined  above  in  sections  3.1.2, 3.1.3,  and  3.1.4,  the  MAESTRO  program  was  used  to  compare 
how  well  the  predicted  values  of  forces  from  the  analytical  model  match  the  wave  balanced 
forces  of  the  MAESTRO  model. 


4.2.3  MAESTRO  Structural  Failure  Modes 


The  specific  failure  modes  analyzed  by  MAESTRO  from  reference  [20]  are  shown  below 
in  Table  15  through  Table  17.  These  failure  modes  are  specific  to  frame  stiffened  structures  such 
as  ships.  The  theoretical  background  behind  these  failure  modes  is  detailed  in  [13]. 


Table  15  -  Panel  Failure  Modes 


PCSF 

Panel  Collapse,  Stiffener  Failure 

PCCB 

Panel  Collapse,  Combined  Buckling 

PCMY 

Panel  Collapse,  Membrane  Yield 

PCSB 

Panel  Collapse,  Stiffener  Flexural/Torsional  Buckling 

PYTF 

Panel  Yield,  Tension  in  Flange 

PYTP 

Panel  Yield,  Tension  in  Plate 

PYCF 

Panel  Yield,  Compression  in  Flange 

PYCP 

Panel  Yield,  Compression  in  Plate 

PSPB  (2  modes) 

Panel  Serviceability,  Plate  Bending 

PFLB 

Panel  Failure  Local  Buckling 

Table  16  -  Frame  Failure  Modes 


FCPH  (3  modes) 

Frame  Collapse,  Plastic  Hinge 

FYCF  (3  modes) 

Frame  Yield,  Compression  in  Flange 

FYTF  (3  modes) 

Frame  Yield,  Tension  in  Flange 

FYCP  (3  modes) 

Frame  Yield,  Compression  in  Plate 

FYTP  (3  modes) 

Frame  Yield,  Tension  in  Plate 

Table  17  -  Girder  Failure  Modes 


GCT 

Girder  Collapse,  Tripping 

GCCF 

Girder  Collapse,  Compression  in  Flange 

GCCP 

Girder  Collapse,  Compression  in  Plate 

GYCF 

Girder  Yield,  Compression  in  Flange 

GYCP 

Girder  Yield,  Compression  in  Plate 

GYTF 

Girder  Yield,  Tension  in  Flange 

GYTP 

Girder  Yield,  Tension  in  Plate 

MAESTRO  measures  the  structural  adequacy  of  each  limiting  failure  mode  using  an 
adequacy  parameter  g(R)  shown  in  (12). 


Instead  of  using  partial  safety  factors  (y)  of  1.25  and  1.5  for  the  serviceability  and 
collapse  failure  modes  respectively,  which  are  consistent  with  ship  structural  design  practice, 
partial  safety  factors  of  1  were  used.  Using  safety  factors  of  one  is  known  as  MAESTRO’s 
Forensic  Mode.  Using  MAESTRO  in  its  Forensic  Mode  allowed  the  analytical  and  MAESTRO 
results  to  be  compared  on  an  equal  basis.  Where  R  is  the  strength  ratio  defined  as  the  fraction  of 
loading  compared  to  the  loading  limit  as  shown  in  (13). 


(13) 


For  each  failure  mode  to  be  satisfied  in  each  loading  condition,  the  adequacy  parameter 
must  be  above  0  to  be  above  the  failure  criteria  for  that  given  failure  mode. 


4.2.4  MAESTRO  Model  Organization 


The  organizational  tree  structure  of  a  general  MAESTRO  Model  is  shown  below  in 
Figure  20. 


Strake  Level 


Figure  20  -  Organizational  Structure  of  a  General  MAESTRO  Model  [20] 


The  overall  MAESTRO  model  is  broken  up  in  to  substructures,  modules,  strakes,  and 
elements.  The  specific  failure  modes  mentioned  above  in  Table  15  through  Table  17  are 
analyzed  at  the  strake  level.  A  strake  consists  of  plates,  stiffeners,  frames,  and  possibly  girders, 
in  which  the  stiffeners,  frames,  and  girders  are  often  T-shapes  for  ship  structures.  This  strake 
evaluation  of  the  failure  modes  is  what  makes  MAESTRO  a  macro-finite  element  program. 

For  the  specific  model  shown  in  Figure  16  through  Figure  18,  the  cross-deck  structure 
between  the  main  and  side  hulls  has  two  substructures  and  five  modules.  While  the  entire 


Figure  21  -  Organization  for  the  Specific  MAESTRO  Model  Analyzed 
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4.2.5  MAESTRO  Verification  of  Analytical  Trimaran  Model  Predicted  Forces 

Once  the  analytical  force  predictions  from  Table  13  and  the  wave  heights  corresponding 
to  5  standard  deviations  above  the  mean  from  Table  14  at  the  applicable  susceptible  wavelengths 
were  obtained,  a  basis  for  comparing  and  testing  the  analytical  model’s  results  to  that  obtained 
by  MAESTRO  in  a  full  quasi-static  linear  theory  wave  balance  was  established.  This 
comparison  was  accomplished  by  applying  the  analytical  model’s  predicted  forces  from  Table  13 
as  point  forces  in  the  finite  element  model  to  the  Stillwater  MAESTRO  model  of  the  trimaran 
shown  in  Figure  16  through  Figure  18,  and  then  measuring  as  an  output  the  failure  modes  shown 
in  Table  15  through  Table  17.  Next,  the  MAESTRO  output  using  the  full  quasi-static  linear 
theory  wave  balance  function  provided  by  the  MAESTRO  program  was  analyzed.  While  the  full 
output  files  from  MAESTRO  are  over  1800  pages,  excerpts  from  these  output  files  are  provided 
in  Appendix  E  through  Appendix  J.  Appendix  E  through  Appendix  J  show  the  adequacy 
parameters  calculated  for  both  the  analytical  force  predictions  and  the  full  MAESTRO  quasi¬ 
static  linear  theory  wave  balance  of  the  specific  strakes  from  Figure  21  for  each  possible  failure 
mode  from  Table  15  through  Table  17  of  the  cross  deck  structure  for  the  trimaran  model .  The 
discussion  of  the  comparison  of  the  analytical  model  to  the  MAESTRO  Model  results  is  included 
in  section  5.2. 


5.1  Statistical  Analysis  of  Analytical  Trimaran  Model  Using  JMP 

The  tabular  results  of  the  analytical  box  shaped  trimaran  model  for  various  main  hull 
lengths,  side  hull  lengths,  side  hull  transverse  spacing,  and  side  hull  longitudinal  placement  from 
section  4.1. 1.7  are  shown  in  Appendix  B.  To  determine  the  relations  between  the  various 
characteristics  of  trimarans  and  the  forces  and  moments  generated  on  a  trimaran’s  cross 
structure,  the  statistical  discovery  program  JMP,  developed  by  the  SAS  institute,  was  used,  JMP 
can  be  used  for  virtually  any  type  of  statistical  analysis.  However,  in  the  case  of  the  analytical 
trimaran  model,  JMP  was  used  to  fit  the  force  and  moment  data  produced  to  three  independent 
variables  in  lengths  and  spacing  to  a  forth  order  polynomial  including  all  the  applicable  cross 
terms  between  the  independent  variables. 

For  each  load  case  the  order  of  polynomial  used  to  fit  the  predicted  data  was  increased 
until  the  R-squared  fit  parameter  and  the  predicted  versus  actual  parameters  were  adequate  to 
accurately  predict  the  data  in  OAppendix  B.  The  entire  fourth  order  fit  equations  and  the  curve 
fitting  statistics  are  shown  in  OAppendix  D.  The  interaction  profiles  between  the  independent 
variables  for  each  case  of  longitudinal  cresting/troughing,  longitudinal  side  hull  twisting,  and 
transverse  sagging/hogging  are  shown  in  the  following  sections.  The  interaction  plots  show  the 
interaction  of  the  independent  variables  between  each  other  and  their  contribution  to  the  overall 
forces.  For  each  polynomial  fit  obtained  an  accompanying  equation  is  stated  to  help  predict  what 
forces  an  early  stage  structural  designer  can  allot  to  the  cross-structure. 


5.1.1  Analytical  Model  Results  for  Longitudinal  Side  Hull  Troughing 


Figure  23  shows  the  side  hull  troughing  forces  as  functions  of  the  three  independent 
variables  of  main  hull  length,  side  hull  length  as  a  fraction  of  main  hull  length,  and  side  hull 
longitudinal  placement  for  amidships  as  a  fraction  of  main  hull  length. 
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Figure  23  -  Longitudinal  Wave  Side  Hull  Troughing  Independant  Variable  Interaction 

The  interaction  plot  shows  that  as  main  hull  length  increases  the  vertical  forces  produced 
due  to  side  hull  troughing  increase.  This  result  is  not  surprising  in  view  of  the  rising  wave  height 
versus  wavelength  characteristic  of  the  statistical  sea  states  described  in  Table  7  and  Table  8, 
since  as  wave  length  of  the  encountered  waves  increases  the  statistically  possible  wave  heights 
increase. 

An  interesting  result  from  the  bottom  row  of  the  interaction  plot  is  that  as  the  longitudinal 
location  of  the  side  hulls  is  moved  aft,  the  trimaran’s  cross  structure  is  not  as  susceptible  to  side 
hull  troughing  as  when  the  side  hull  is  placed  in  the  longitudinally  symmetric  amidships  position. 


Another  result  from  the  middle  row  of  Figure  23  is  that  vertical  forces  due  to  longitudinal 
troughing  are  only  mildly  sensitive  to  changes  in  the  length  of  the  side  hull,  but  lessen  as  the 
length  of  the  side  hull  increases.  Finally,  the  interaction  between  variables  can  be  studied  as 
well.  For  example,  having  a  smaller  side  hull  length  produces  higher  vertical  forces  when  the 
side  hulls  are  placed  further  aft  than  when  the  same  side  hull  is  placed  amidships. 

One  final  result  from  the  side  hull  troughing  data  is  the  fit  equation  provided  in  (14). 

This  equation  gives  the  downward  force  as  a  function  of  three  independent  variables.  For  clarity 
the  equation  (14)  is  only  shown  to  second  order  with  the  full  fourth  order  equation  included  in 
Appendix  D. 
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5.1.2  Analytical  Model  Results  for  Longitudinal  Side  Hull  Cresting 

The  interaction  profiles  for  the  upward  vertical  forces  produced  during  side  hull  cresting 
are  shown  below  in  Figure  24. 
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Figure  24  -  Longitudinal  Wave  Side  Hull  Cresting  Independant  Variable  Interaction 

The  interaction  plots  for  side  hull  cresting  is  analogous  to  that  of  side  hull  troughing,  and 
the  relationships  between  the  independent  variables  are  essentially  the  same  as  that  of  side  hull 
troughing. 
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5.1.3  Analytical  Model  Results  for  Longitudinal  Positive  Phase  Twisting 

The  polynomial  equation  fit  results  for  the  positive  phase  twisting  of  the  cross  structure  in 
longitudinal  waves  is  shown  below  in  the  interaction  plot  of  Figure  25. 
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Figure  25  -  Longitudinal  Wave  Side  Hull  Positive  Twisting  Independant  Variable  Interaction 


Like  the  interaction  plots  for  longitudinal  cresting  and  troughing,  the  moments  produced 
in  positive  phase  twisting  increases  markedly  as  main  hull  length  increases  as  is  expected  from 
the  wavelengths  of  the  susceptible  sea  states.  From  the  bottom  row  of  Figure  25  it  is  observed 
that  the  longitudinal  placement  of  the  side  hulls  has  virtually  no  effect  on  the  positive  phase 
twisting  moments.  However,  from  the  second  row  of  the  interaction  plot,  it  can  be  seen  that  as 
the  length  of  the  side  hull  increases  the  twisting  moment  increases  drastically.  This  result  is  not 
unreasonable  since  the  length  of  the  side  hull  increases  the  length  of  the  moment  arm  for  positive 
phase  twisting. 


The  partial  fit  equation  shown  for  positive  phase  twisting  is  shown  equation  (16)  with  the 
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5.1.4  Analytical  Model  Results  for  Longitudinal  Negative  Twisting 

The  polynomial  equation  fit  results  for  the  negative  phase  twisting  of  the  cross  structure 
in  longitudinal  waves  is  shown  below  in  the  interaction  plot  of  Figure  26. 
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Figure  26  -  Longitudinal  Wave  Side  Hull  Negative  Twisting  Independant  Variable  Interaction 

The  results  for  negative  phase  twisting  are  analogous  to  positive  phase  twisting 
maintaining  similar  trends  between  the  independent  variables. 
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The  partial  fit  equation  shown  for  positive  phase  twisting  is  shown  equation  (17)  with  the 
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5.1.5  Analytical  Model  Results  for  Transverse  Sagging 

The  results  for  transverse  sagging  are  shown  for  the  interaction  plot  below  in  Figure  27. 
For  clarity  the  values  in  Appendix  D  were  converted  to  utilize  transverse  dimensions  stated  in 
main  and  side  hull  beams  from  the  fixed  relations  in  Table  11.  The  interaction  plots  indicate 
that,  as  suspected,  the  upward  forces  on  the  cross  structure  due  to  a  transverse  sagging  condition 
increases  as  the  outer  hull  spacing  increases.  Another  result  from  the  interaction  profile  is  that  as 
the  side  hulls  beam  widens  then  forces  on  the  cross  structure  lessens.  Finally,  as  experienced  in 
the  previous  loading  cases  as  the  size  of  the  overall  ship  increases  the  forces  experienced 
increase  due  to  the  statistical  nature  of  the  sea  states. 
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Figure  27  -  Transverse  Sagging  Wave  Side  Hull  Independant  Variable  Interaction 


The  partial  fit  equation  shown  for  transverse  sagging  is  shown  equation  (1 8)  with  the  full 
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5.1.6  Analytical  Model  Results  for  Transverse  Hogging 

A  polynomial  fit  equation  for  the  downward  forces  produced  during  transverse  hogging 
was  unable  to  be  fitted  to  the  data  in  Appendix  B  even  using  forth  order  terms.  This  situation 
arises  from  the  fact  that  there  is  a  discontinuity  in  the  downward  force  applied  in  a  transverse 
hogging  condition.  When  the  transverse  wave  height  is  small  enough  that  the  outer  hulls  do  not 
leave  the  water  the  loss  of  buoyancy  due  to  wave  height  is  linear.  However,  when  the  wave 
height  becomes  large  enough,  the  side  hulls  leave  the  water  and  the  loss  of  buoyancy  after  that 
point  regardless  of  the  wave  height  remains  constant.  Hence,  a  good  fit  for  an  overall  downward 
force  due  to  transverse  hogging  over  a  large  range  of  values  was  not  able  to  be  found. 
Alternatively,  the  structural  designer  of  a  trimaran  should  design  to  the  standard  of  the  side  hulls 
completely  broaching  the  water  as  is  discussed  in  references  [5]  through  [10]. 


5.2  Comparison  and  Discussion  of  Analytical  and  MAESTRO  Results 

An  example  of  the  results  of  the  MAESTRO  structural  analyses  of  the  analytical  force 
predictions  compared  to  the  full  quasi-static  linear  wave  theory  balance  analysis  from  4.2.5  is 
shown  below  in  Table  18  and  Table  19  with  the  full  results  for  all  load  cases  and  modules 
included  in  Appendix  E  through  Appendix  J. 


Table  18  -  Adequacy  Parameters  with  analytically  Predicted  Forces  (Longitudinal  Troughing) 

INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  1  OF  SUBSTR.  1 


STRAKE 

PCSF 

PCCB 

PCMY 

PCSB 

PYTF 

PYTP 

PYCF 

PYCP 

PSPBT 

PSPBL 

PFLB 

1 

0.826 

0.927 

0.965 

0.965 

1.000 

• 

1.000 

2 

0.883 

0.978 

0.939 

0.945 

1.000 

1.000 

0.882 

3 

0.971 

0.986 

1.000 

1.000 

0.923 

4 

0.939 

0.955 

0.991 

0.995 

1.000 

1.000 

0.932 

5 

0.979 

0.983 

1.000 

1.000 

0.979 

6 

0.957 

0.976 

0.996 

0.996 

1.000 

1.000 

0.901 

7 

0.913 

1.000 

1.000 

0.906 

+ - - - - - - - + 


Table  19  -  Adequacy  Parameters  with  MAESTRO  Wave  Balance  (Longitudinal  Troughing) 

INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  1  OF  SUBSTR.  1 


|  STRAKE 

PCSF 

PCCB 

PCMY 

PCSB 

PYTF 

PYTP 

PYCF 

PSPBT 

PSPBL  i 

PFLB 

i  1 

0.848 

0.943 

0.937 

0.939 

1.000 

1.000 

0.969 

0.969 

1.000 

1.000 

0.831 

1  2 

0.898 

0.980 

0.944 

0.918 

1.000 

1.000  | 

0.948 

0.948 

1.000 

1.000 

0.897 

j  3 

0.972 

0.997 

0.988 

0.986 

1.000 

1.000 

0.991 

0.991 

1.000 

1.000 

0.940 

i  4 

0.949 

0.987 

0.964 

0.991 

1.000 

1.000 

0.995 

0.995 

1.000 

1.000 

0.944 

1  5 

0.977 

1.000 

0.983 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.977 

I  6 

0.967 

0.986 

0.983 

0.995 

1.000 

1.000 

0.997 

0.997 

1.000 

1.000 

0.923 

|  7 

+ - 

0.918 

0.967 

0.945 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.908 

Appendix  E  through  Appendix  J  show  the  adequacy  parameters  for  the  panel  and  frame 
failure  modes  for  each  strake  in  the  cross-deck  structure  from  Figure  21.  From  the  way  that  the 
adequacy  parameter  is  defined,  a  lower  adequacy  parameter  means  that  the  evaluated  strake  is 
closer  to  failure  for  the  analyzed  failure  mode. 

Since  a  lower  adequacy  parameter  means  closer  to  failure,  for  the  analytical  model  to  be  a 
conservative  estimate  of  the  cross-deck  structural  loading  of  the  trimaran,  every  adequacy 
parameter  for  the  analytical  model  must  be  less  than  the  corresponding  MAESTRO  model 
adequacy  parameter.  While  the  adequacy  parameters  in  Appendix  E  through  Appendix  J  do 


compare  to  a  first  order,  the  analytical  model  is  not  consistently  the  conservative  estimate  for 
every  failure  mode  for  every  strake  of  the  cross-deck  structure. 

The  fact  that  the  analytical  model  does  not  conservatively  predict  the  failure  adequacy 
parameters  is  not  surprising.  Many  simplifications  to  the  analytical  model  were  made  in  order  to 
make  it  possible  to  investigate  the  basic  nature  of  global  loading  of  the  cross-deck  structural 
forces  experienced  by  trimarans  over  a  large  range  of  hull  sizes  and  hull  spacing  configurations. 
The  simplification  most  affecting  the  performance  of  the  analytical  model  compared  to  the  fully 
wave  balanced  MAESTRO  model  was  the  assumption  that  the  hulls  were  rigid  and  rigidly 
connected.  The  hull  deflections  in  the  realistic  fully  wave  balanced  models  were  the  primary 
source  of  the  difference  between  the  two  models. 

While  it  is  unfortunate  that  the  analytical  model  is  not  a  truly  conservative  estimate  of  the 
cross-deck  structural  loading  of  trimarans  for  every  loading  case  and  failure  mode,  the  results  of 
the  two  methods  did  compare  to  a  first  order  of  magnitude  with  each  other.  Since  the  analytical 
model  can  predict  forces  to  a  first  order  of  magnitude,  it  can  be  a  useful  tool  for  early  stage 
design  estimates  of  trimaran  cross  deck  structural  loading  as  functions  of  the  basic  hull 
parameters  and  dimensions  restated  again  in  Table  20  and  Table  21. 


Table  20  -  Structural  Loading  from  Sections  3.1.2, 3.1.3,  and  3.1.4 


Jgto  _ 

the  vertical  force  on  cross-structure  in  longitudinal  waves  (troughing  and  cresting) 

MBig 

the  moment  on  cross-structure  in  longitudinal  waves  (longitudinal  twisting) 

FB,,. 

the  vertical  force  in  transverse  waves  (transverse  hogging/sagging) 

Table  21  -  Relevant  Trimaran  Parameters  Affecting  Design  Loading 


Bm 

the  main  hull’s  beam 

Lm 

the  main  hull’s  length 

Tm 

the  main  hull’s  draft 

Fm 

the  main  hull’s  freeboard  (also  side  hull  freeboard) 

Bs 

the  side  hull’s  beam 

u 

the  side  hull’s  length 

Ts 

the  side  hull’s  draft 

X, 

the  longitudinal  position  of  the  side  hull  wrt  main  hull  amidships 

Ys 

the  transverse  position  of  the  side  hull  wrt  main  hull  centerline 

Since  the  goal  of  this  work  was  to  provide  the  trimaran  designer  with  a  way  to  estimate 
the  cross-deck  structural  loading,  this  work  is  considered  a  qualified  success. 


Chapter  6  Conclusions 


This  main  product  of  this  work  was  the  curve  fits  in  Appendix  D  that  predict  trimaran 
cross-deck  structure  loading  in  applicable  load  cases  of  longitudinal  troughing/cresting, 
longitudinal  positive/negative  twisting,  and  transverse  hogging/sagging.  These  fitted  curve  of 
design  loadings  allows  an  initial  design  stage  loading  estimate  for  cross  deck  structural  loading 
given  general  characteristics  of  length  and  spacing  of  a  trimaran’s  hulls.  The  actual  equations 
derived  from  the  analytical  model  are  useful  for  first  approximations  of  loading  to  the  cross-deck 
structure  of  a  trimaran  but  are  not  necessarily  always  conservative.  Flexure  of  the  main  hull  of 
the  trimaran  is  the  main  cause  of  the  analytical  model’s  un-conservative  structural  loading 
predictions. 

A  concurrent  result  of  the  fourth  order  polynomial  fitted  equations  were  the  interaction 
profile  plots  shown  in  Figure  23  through  Figure  27.  The  interaction  profile  plots  show  how  the 
main  design  variables  of  the  trimaran  interact  with  each  other  to  affect  the  cross-deck  global 
structural  loading.  These  plots  are  a  useful  visual  qualitative  tool  to  determine  which  trimaran 
configurations  experience  less  cross-deck  structural  loading. 

The  fourth  order  polynomial  curve  fitted  equations  and  the  interaction  profile  plots 
combined  with  other  characteristics  of  good  trimaran  design  including  stability,  roll,  and 
resistance  characteristics  will  aid  the  trimaran  ship  designer  in  optimizing  an  overall  trimaran 
ship  design. 


Chapter  7  Recommendations  for  Future  Work 


The  first  and  probably  most  important  area  to  continue  further  work  in  the  area  of 
trimaran  cross-deck  structural  design  would  be  to  have  an  analytical  model  that  accounted  for 
forward  speed  effects  of  the  ship  moving  through  the  water.  While  this  is  not  extremely 
complicated  to  perform,  it  was  not  included  in  this  work  simply  because  MAESTRO  does  not 
have  the  capability  verify  the  results,  and  the  entire  goal  of  this  work  was  to  derive  an  analytical 
model  for  cross-deck  structural  loading  that  could  be  compared  with  a  more  rigorous  analysis 
such  as  finite  element  analysis. 

An  improvement  to  the  analytical  model  in  this  work  would  be  to  account  for  side  hull 
flare,  slamming  of  waves  into  the  cross  deck  structure,  and  incident  waves  encountered  at 
oblique  angles.  Accounting  for  flare  in  the  side  hulls  would  be  beneficial  because  both  [1]  and 
[19]  reference  the  need  to  have  side  hull  flare  for  stability  and  sea-keeping  reasons.  Using 
triangular  side  hulls  instead  of  box  barge  side  hulls  would  be  a  relatively  simple  correction  to  the 
analytical  model  and  could  be  done  to  make  side  hull  flare  angle  an  additional  parameter  that 
could  be  varied  while  calculating  design  forces.  Wave  slamming  of  cross-deck  structures  and 
obliquely  angled  waves  encounters  on  the  other  hand  is  more  complicated  and  would  require 
extensive  modifications  to  the  current  analytical  model. 

Finally,  general  recommendations  for  study  in  the  area  of  trimaran  structural  design  in 
general  would  be  to  determine  the  dynamic  whipping  response  of  the  relatively  long  and  slender 
hulls  that  are  characteristic  of  trimarans.  An  accompanying  topic  to  the  whipping  response  of 
slender  hulls  would  be  an  investigation  into  active  structural  control  to  include  the  cost  and 
power  requirements  for  very  high  speed  applications. 
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Appendix  A.  MathCAD  Analytical  Model 


All  the  equations  and  formula  required  to  complete  the  analytical  model  of  the  “box”  trimaran 
are  included  in  this  Appendix.  The  essential  explanations  of  the  working  of  the  analytical  model 
are  included  in  section  4.1.1.  Where  appropriate,  each  applicable  section  refers  to  the  main 
text’s  explanation. 

Explained  in  section  4. 1.1.1 
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Explained  in  section  4. 1.1. 2 

B:=(rows(ssorig)-l) 

D:=(cols(ssorig)-l) 

t  :=  1..2B| 

r  :=  1..D 

:=  submatrix^ss  orjg> 0,0,0, 19)T-sec 

hwt:=  (l  ~  -5)  m 

(This  equation  explained  in  Section  4. 1.1. 3) 
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Explained  in  Section  4,1. 1.3 
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motion_tr|p ,  q  ,<j>,  Y  s  ,T 


Explained  in  Section  4. 1.1.5 
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design_lg_hog_FB  s(xs>xcf>Ts,L.s,Bs>Fm,Tm,Bm,Lin) : 
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design_twist_MB_pos  (X$,X cf ,T S,LS ,B s ,F m,T m,B m,L^ := 
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design  4~  0 

q  4—  q_twist_MB_pos  (xg,Xcf/rs,Ls,Bs,Fm,Tm,Bm,L  J 

pmaXq 

sum  4—  >  ss 

i-t  p.q 

P  =  1 

for  p  s  l..pmax 


if  ss  *  0 

p.q 


F  4—  MB 


lg^p’!,’^lg0’Xs’Xcf>Ts*Ls>Bs,Fm’Tm’Bm’I'in] 


p.q 


mean  4-  mean  +  |f|  — — 
sum 


m2  4—  m2+  (F) 


„  ss 
2  p.q 


design  4—  mean 
design 


+  m2  - 


(mean) 


q_twist_MB_neg  (X s, X cf,T s ,LS ,B s ,F m,T m, B m,L  J  := 


‘l  max  4-0 
F  4-  0 

for  q 6  l..dmax_lg(Ln^ 
p  4-  pmax^ 

Ftemp  4-  iMBig^P.q.^ig^.Xj.X^.Tj.Lj.Bj.Fjjj.T^Bnj.Ljjjjl 
if  Ftemp  >  F 
hmax<-‘l 
I F  4—  Ftemp 


design_twist_MB_neg  (Xs,Xcf,Ts,Ls,Bs>Fm,T := 


mean  4-  0 
m2  4-  0 
design  4-  0 

q  4—  q_twist_MB_neg  (xs’Xcf>Ts’Ls,Bs’pm,Tm>Bm'Ln| 

pmaXq 

sum  4—  >  ss 

i-t  P.q 

p  =  l 

for  p  e  l„pmax 


if  ss  *■ 0 

p.q 


F  <—  MB 


Ig^P.q.^lg^.Xg.Xcf.Tj.Lj.Bj.Fjjj.TjjpBnj.lL,,^ 


mean  4— meant  |f|  - 


p.q 


m2  4—  m2  +  (F) 


„  ss 

2  p.q 


design  4 — [mean  +  (mean)^] 

design 


m 
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dmax_ti^Ysj  := 


dmax_tr  4—  0 
for  n  g  1..D 

dmax_tr  4—  n  if  Xn  <  2-Yg 
dmax  tr  4—  dmax  tr  +  3 


q_tr_sag_FB  (' Y  , ,  T  s ,  L  g ,  B  g ,  F  m,  T  m,  B  m,  L  J  := 


max' 

F  4—  0 

for  qe  1.. dmax_tr^Ys) 


p  4—  pmax 


Ftemp  4— 


FB 


if  Ftemp  >  F 
^  max*-  9 
F  <—  Ftemp 


1  max 


design_tr_sag_FB  s(Ys,Ts,Ls>Bs,Fm,Tm,Bm,Lln)  := 


mean  4—  0 
m2  4-  0 
design  4—  0 

q  q_tr_sag_FB  (Ys,Ts,Ls,Bs,Fm,Tm,Bm,LJIJ 

pmaxq 

sum  4—  >  ss 

p.q 

p  =  1 

for  pel.,  pmax 


if  ss 

p.q 


F  4—  FB 


mean  <—  mean 


+  |F|- 


p.q 


m2  4—  m2  +  (F) 


„  ss 

2  p.q 


design  4—  [mean  +  (3->/m2-  (mean)'*] 
design 


79 


q_tr_hog_FB  (Ys,Ts,Ls,Bs.Fm»Tm,Bm.L  J . 


^  max  ®  ' 

F  <—  0 

for  qe  l..dmax_ti|Ys| 
p  <—  pmax^ 

Ftemp  <—  jFBtr^p,q,^>tr^,Ys,Ts,Ls,Bs,Fin,Tln,Bm>Lin’j| 

if  Ftemp  >  F 

Q  max  Q 
F  <—  Ftemp 

I  max 


m 


desigti_tr_hog_FB  ,(Y8.T#,L5,Bs,Fm,Tffi.Bm.L  J  := 


mean  <—  0 
m2 1 —  0 
design  <—  0 

Q  <—  q_tr_hog_FB  (Yj.Tj.Lj.Bj.Fjjj.Tjjj.Bjjj.Ljj^ 

pmaxq 

sum  <—  >  ss 

£->  p.q 

p  =  1 

for  pe  l..pmax 


if  ss  *  0 

p.q 


F  <—  FB 


t/p>ci.<t,tr, 
V  c 


’Ys>Ys.Ls.Bs,Fm,Tm,Bm,Lm , 


mean  <—  mean 


+  |f|- 


ss 


p.q 


sum 


m2  <—  m2  +  (F) 


„  ss 
2  p.q 

sum 


design  < — [mean  +  ji-J m2  -  (mean)' 
design 


Explained  in  Section  4. 1.1. 7 

file_root  :=  "C:\Documents  and  Settings\jlrhoads\My  DocumentsNrhoads  thesis  data\" 
fiIe_base(L_)  :=num2str — 

'  v  V  m  J  file_ext:=  ”.prn"  cc(a,b)  :=concat(a,b)  | 

fullname Jgf L  :=  File_root  cc  j^fileJjase^L^)  cc  ("_lg"  cc  file_ex()JJ  | 

fulIname_tr(L  ^  :=  fde_root  cc  j^file^ase^L^  cc  ("_tr"  cc  file_ex§JJ  | 
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Bn 


for  j  e  0..  12 
|Lm(— (60+20 j)m 

14 

2 

•Dm-’ 

data  jg4—  (0  000000  0) 
data  tr  «—  (0  0  0  0  0  0) 
WRITEPRI^fulIname Jg  (L  J ,  data 
WRITEPRr^ful!naine_tr  (L  ^ ,  data 
for  ne  0..5 

L_  4—  Lm-(0.25+  n-.05) 


LmBn 

",  Ls 


Bs+-Ts 
for  i  e  0..  5 


"(t-t) 


xs  xmc'' 

x  BmLm‘B‘m+  ^Bs*Ls'Xs 
®m^m+  ^Bs'Ls 
Lm 


data  i. 


0,0  m 


data .  4 - 

«0.1  m 

,  Bs 

Jg0,2  m 

„  X* 

data  ip  4 - 

^0,3  m 

design Jg.sag.FB  s(xs,Xcf,Ts,Ls,Bs,Fm,Tm,Bm,Lnl) 

data  .  4 - — 

Ig0,4  N 

design Jg_hog_FB  s(X!,Xcf,Ts,Ls,Bs,Fm,TnrBm,L  j 

data  k  4 - — 

^0,5  N 

design_twist_MB_pos  (xs,Xcf,Ts,Ls,Bs,Fm,Tm,Bm,Lm) 

data  4 - — - 

*"0,6  N-m 

design_twist.MB.neg  (xs,Xcf,Ts,Ls,Bs,Fm,TIn,Bm,Lm) 

data  k  4 - - 

^0,7  N-m 

APPENDPRl^fullnameJg  (L  ^ ,  data  jg) 
for  ie  0..5 

Ys<-Bm-(l  +  iO.I) 


0,0  m 

L_ 


data 


trc 


0, 1  m 
^s 

‘0,2  m 

Yc 


tin 


0,3  m 

design.tr_sag.FB  s(Ys,Ts,Ls,Bs,Fm,Tm,Bm,L J 

data  fr  4 - - - - 

**0,4  N 

design_tr_hog_FB  s(Ys,Ts,Ls,Bs,Fm,Tm,Bm,L J 

data  fr  4-  — — - 

ir0,5  N 

APPENDPRI^fuIInarne_tr  (l  ^ ,  data 


"done" 
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Appendix  B.  MathCAD  Analytical  Results  Tables 
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5.71  E+05 

80 

32 

1.278 

8 

6.34E+05 

80 

32 

1.278 

8.571 

6.93E+05 

80 

36 

1.205 

5.714 

5.22E+05 

80 

36 

1.205 

6.286 

5.59E+05 

80 

36 

1.205 

6.857 

5.76E+05 

80 

36 

1.205 

7.429 

6.00E+05 

80 

36 

1.205 

8 

6.66E+05 

80 

36 

1.205 

8.571 

7.28E+05 

80 

40 

1.143 

5.714 

5.46E+05 

80 

40 

1.143 

6.286 

5.85E+05 

80 

40 

1.143 

6.857 

6.02E+05 

80 

40 

1.143 

7.429 

6.27E+05 

80 

40 

1,143 

8 

6.96E+05 

80 

40 

1.143 

8.571 

7.61  E+05 

F_tr_hog  (N) 
-3.98E+05 
-4.21  E+05 
-4.32E+05 
-4.29E+05 
-4.13E+05 
-3.83E+05 
-4.26E+05 
-4.52E+05 
-4.63E+05 
-4.60E+05 
-4.43E+05 
-4.12E+05 
-4.51  E+05 
-4.79E+05 
-4.92  E+05 
-4.88E+05 
-4.69E+05 
-4.36E+05 
-4.75E+05 
-5.04E+05 
-5.18E+05 
-5.14E+05 
-4.94E+05 
-4.58E+05 
-4.96E+05 
-5.28E+05 
-5.42E+05 
-5.38E+05 
-5.16E+05 
-4.79E+05 
-5.17E+05 
-5.49E+05 
-5.65E+05 
-5.61  E+05 
-5.38E+05 
-4.98E+05 


100 

25 

2.02 

0 

-1.14E+06 

1.17E+06 

5.59E+06 

-5.59E+06 

100 

25 

2.02 

-7.5 

-1.07E+06 

1 .09E+06 

5.43E+06 

-4.73E+06 

100 

25 

2.02 

-15 

-8.95E+05 

8.99E+05 

4.82E+06 

-4.85E+06 

100 

25 

2.02 

-22.5 

-6.95E+05 

6.95E+05 

4.70E+06 

-4.87E+06 

100 

25 

2.02 

-30 

-5.74E+05 

5.74E+05 

4.63E+06 

-5.05E+06 

100 

25 

2.02 

-37.5 

-5.80E+05 

5.80E+05 

4.70E+06 

-5.52E+06 

100 

30 

1.844 

0 

-1.15E+06 

1.18E+06 

8.03E+06 

-8.03E+06 

100 

30 

1.844 

-7 

-1.09E+06 

1.11E+06 

7.83E+06 

-7.80E+06 

100 

30 

1.844 

-14 

-9.14E+05 

9. 22 E+05 

7.15E+06 

-6.16E+06 

100 

30 

1.844 

-21 

-7.10E+05 

7.12E+05 

7.47E+06 

-6.23E+06 

100 

30 

1.844 

-28 

-5.65E+05 

5.66E+05 

6.55E+06 

-7.44E+06 

100 

30 

1.844 

-35 

-5.34E+05 

5.34E+05 

6.76E+06 

-8.24E+06 

100 

35 

1.707 

0 

-1.12E+06 

1.16E+06 

1.05E+07 

-1.05E+07 

100 

35 

1.707 

-6.5 

-1.07E+06 

1.09E+06 

1.03E+07 

-1.02E+07 

100 

35 

1.707 

-13 

-9.06E+05 

9.20E+05 

9.77E+06 

-9.51  E+06 

100 

35 

1.707 

-19.5 

-7.07  E+05 

7.1 1  E+05 

1.02E+07 

-6.99E+06 

100 

35 

1.707 

CD 

CM 

1 

-5.45E+05 

5.46E+05 

8.47E+06 

-1 .01  E+07 

100 

35 

1.707 

-32.5 

-4.79E+05 

4.79E+05 

8.58E+06 

-1.13E+07 

100 

40 

1.597 

0 

-1.19E+06 

1.11E+06 

1.29E+07 

-1 .29  E+07 

100 

40 

1 .597 

-6 

-1 .02E+06 

1.05E+06 

1.26E+07 

-1.25E+07 

100 

40 

1.597 

-12 

-8.75E+05 

8.95E+05 

1.25E+07 

-1.17E+07 

100 

40 

1.597 

-18 

-6.85E+05 

6. 94 E+05 

1.31  E+07 

-1.25E+07 

100 

40 

1.597 

-24 

-5.15E+05 

5.18E+05 

1.35E+07 

-1.27E+07 

100 

40 

1.597 

-30 

-4.18E+05 

3.21  E+05 

1.01  E+07 

-1.43E+07 

100 

45 

1.506 

0 

-1.15E+06 

1.05E+06 

1.49E+07 

-1.49E+07 

100 

45 

1.506 

-5.5 

-1 .09E+06 

1.13E+06 

1.46  E+07 

-1.45E+07 

100 

45 

1.506 

-11 

-8.20E+05 

9.51  E+05 

1.52E+07 

-1.35E+07 

100 

45 

1.506 

-16.5 

-6.49E+05 

7.22E+05 

1.59E+07 

-1 .50E+07 

100 

45 

1.506 

-22 

-4.77E+05 

4.75E+05 

1.65E+07 

-1.52E+07 

100 

45 

1.506 

-27.5 

-3.76  E+04 

3.03E+05 

1 .68E+07 

-1.72E+07 

100 

50 

1.429 

0 

-1.08E+06 

9.57E+05 

1 .63E+07 

-1.63E+07 

100 

50 

1.429 

-5 

-1.03E+06 

1.08E+06 

1 .70  E+07 

-1.59E+07 

100 

50 

1.429 

-10 

-2.19E+05 

9.27E+05 

1 .76  E+07 

-1.47E+07 

100 

50 

1.429 

-15 

-1 .30E+05 

6.79E+05 

1 .84  E+07 

-1.72E+07 

100 

50 

1.429 

o 

CM 

* 

-7.48E+04 

4.88E+05 

1.91  E+07 

-1.73E+07 

100 

50 

1.429 

-25 

-7.57E+04 

3.02E+05 

1 .97E+07 

-1.96E+07 

nsn 

■5531 

■HfiSl 

KI3G31 

BSBYQjl 

100 

25 

2.02 

7.143 

6.58E+05 

-6.41  E+05 

100 

25 

2.02 

7.857 

7.60E+05 

25 

2.02 

8.571 

8.54E+05 

k mm 

KB 

25 

2.02 

9.286 

9.38E+05 

mm mm 

100 

25 

2.02 

10 

1.01  E+06 

-9.99E+05 

100 

25 

2.02 

10.71 

1.06E+06 

-1.05E+06 

100 

30 

1.844 

7.143 

7.13E+05 

-6.96E+05 

100 

30 

1.844 

7.857 

8.23E+05 

-6.67E+05 

mm 

30 

1.844 

8.571 

9.25E+05 

-9.20E+05 

— 

30 

1.844 

9.286 

1.02E+06 

-1.01  E+06 

100 

30 

1.844 

10 

1.09  E+06 

-1.08E+06 

100 

30 

1.844 

10.71 

1.15E+06 

-1.13E+06 

100 

35 

1.707 

7.143 

7.62  E+05 

-7.44E+05 

100 

35 

1.707 

7.857 

8.80E+05 

-7.14E+05 

KB 

35 

1.707 

8.571 

9.89E+05 

-9.81  E+05 

100 

35 

1.707 

9.286 

1.09 E+06 

-1.07E+06 

100 

35 

1.707 

10 

1.17E+06 

-1.14E+06 

100 

35 

1.707 

10.71 

1.23E+06 

-1.19E+06 

100 

40 

1.597 

7.143 

8.07E+05 

-7.81  E+05 

100 

40 

1.597 

7.857 

9.32E+05 

-7.53E+05 

100 

40 

1.597 

8.571 

1 .05  E+06 

-6.79E+05 

100 

40 

1.597 

9.286 

1.15E+06 

-1.12E+06 

100 

40 

1.597 

10 

1.24E+06 

-1.19E+06 

100 

40 

1.597 

10.71 

1.30E+06 

-1.24E+06 

100 

.  45 

1.506 

7.143 

8.48E+05 

-8.14E+05 

100 

45 

1.506 

7.857 

9.80E+05 

-7.86E+05 

100 

45 

1.506 

8.571 

1.10E+06 

-7.1 2  E+05 

100 

45 

1.506 

9.286 

1.21  E+06 

-1.17E+06 

100 

45 

1.506 

10 

1.30E+06 

-1 .24  E+06 

100 

45 

1.506 

10.71 

1.37E+06 

-1.29E+06 

KB 

50 

1.429 

7.143 

8.87E+05 

-8.46E+05 

KB 

50 

1.429 

7.857 

1.02E+06 

-8.17E+05 

100 

50 

1.429 

8.571 

1.15E+06 

-7.42  E+05 

100 

50 

1.429 

9.286 

1 .26E+06 

-1.21  E+06 

100 

50 

1.429 

10 

1.36E+06 

-1 .28E+06 

100 

50 

1.429 

10.71 

1.43E+06 

-1.34E+06 

lOQjl 

Bs  (m) 

ilviniiB 

rnZMmm 

'mMmfim  - '1 

120 

30 

2.424 

0 

-1.56E+06 

1.94E+06 

1.03E+07 

-1.03E+07 

120 

30 

2.424 

-9 

-1 .50E+06 

1 .51 E+06 

1.01  E+07 

-9.99E+06 

120 

30 

2.424 

-18 

-1.45E+06 

1.35E+06 

9.79E+06 

-9.30E+06 

120 

30 

2.424 

-27 

-1.11E+06 

1.11  E+06 

9.43E+06 

-8.17E+06 

120 

30 

2.424 

-36 

-8.82E+05 

8.83E+05 

9.24E+06 

-9.55E+06 

120 

30 

2.424 

-45 

-9.39E+05 

9.38E+05 

9.84E+06 

-1.06E+07 

120 

36 

2.213 

0 

-1.95E+06 

1.98E+06 

1.50E+07 

-1 .50E+07 

120 

36 

2.213 

-8.4 

-1.83E+06 

1.85E+06 

1.49E+07 

-1.38E+07 

120 

36 

2.213 

-16.8 

-1.51E+06 

1 .33E+06 

1.44E+07 

-1 .28E+07 

120 

36 

2.213 

-25.2 

-1.13E+06 

1.10E+06 

1.39E+07 

-1 .21  E+07 

120 

36 

2.213 

-33.6 

-8.75E+05 

8.75E+05 

1.35E+07 

-1 .41  E+07 

120 

36 

2.213 

-42 

-8.53E+05 

8.53E+05 

1 .34E+07 

-1 .58E+07 

120 

42 

2.049 

0 

-1.92E+06 

1.97  E+06 

2.01  E+07 

-2.01  E+07 

120 

42 

2.049 

-7.8 

-1.81  E+06 

1.85E+06 

2.00E+07 

-1.72E+07 

120 

42 

2.049 

-15.6 

-1.52E+06 

1.53E+06 

1.94E+07 

-1.58E+07 

120 

42 

2.049 

-23.4 

-1.15E+06 

1.15E+06 

1 .87E+07 

-1 .48E+07 

120 

42 

2.049 

-31.2 

-8.56E+05 

8.57E+05 

1.81  E+07 

-1.91  E+07 

120 

42 

2.049 

-39 

-7.59E+05 

7.59E+05 

1 .79E+07 

-2.16E+07 

120 

48 

1.917 

0 

-1 .86E+06 

1.92E+06 

2.53E+07 

-2.53E+07 

120 

48 

1.917 

CM 

h-i 

1 

-1 .76E+06 

1.81  E+06 

2.51  E+07 

-2.48E+07 

120 

48 

1.917 

-14.4 

-1.49E+06 

1.52E+06 

2.45E+07 

-2.40E+07 

120 

48 

1.917 

-21.6 

-1.14E+06 

1.15E+06 

2.36E+07 

-2.30E+07 

120 

48 

1.917 

-28.8 

-8.28E+05 

8.30E+05 

2.28E+07 

-2.60E+07 

120 

48 

1.917 

-36 

-6.66E+05 

6.67E+05 

2.23E+07 

-2.75E+07 

120 

54 

1.807 

0 

-1 .88E+06 

1.83E+06 

3.00E+07 

-3.00E+07 

120 

54 

1.807 

CD 

CD 

1 

-1 .78E+06 

1.73E+06 

2.99E+07 

-2.95E+07 

120 

54 

1.807 

-13.2 

-1 .43  E+06 

1.47E+06 

2.92E+07 

-2.85E+07 

120 

54 

1.807 

-19.8 

-1.11  E+06 

1.12E+06 

2.82E+07 

-2.73E+07 

120 

54 

1.807 

-26.4 

-7.91  E+05 

7.97E+05 

2.71  E+07 

-3.12E+07 

120 

54 

1.807 

-33 

-9.56E+04 

5.82E+05 

2.64E+07 

-3.30E+07 

120 

60 

1.714 

0 

-1.78E+06 

1 .70E+06 

3.39E+07 

-3.39E+07 

120 

60 

1.714 

-6 

-1.70E+06 

1 .62E+06 

3.38E+07 

-3.33E+07 

120 

60 

1.714 

-12 

-1.34E+06 

1 .51  E+06 

3.31  E+07 

-3.22E+07 

120 

60 

1.714 

-18 

-3.60E+05 

1.16E+06 

3.53E+07 

-3.09E+07 

120 

60 

1.714 

-24 

-2.12E+05 

7.87E+05 

3.79E+07 

-3.58E+07 

120 

60 

1.714 

o 

CO 

1 

-7.88E+04 

4.68E+05 

4.00E+07 

-3.78E+07 

mm 

EX&3] 

WSSSi 

120 

30 

2.424 

8.571 

1.18E+06 

-1.18E+06 

120 

30 

2.424 

9.429 

1.32E+06 

-1.32E+06 

120 

30 

2.424 

10.29 

1.44E+06 

-1.43  E+06 

120 

30 

2.424 

11.14 

1 .52E+06 

-1.51  E+06 

120 

30 

2.424 

12 

1.56E+06 

-1.55E+06 

120 

30 

2.424 

12.86 

1.56E+06 

-1.55E+06 

120 

36 

2.213 

8.571 

1.28E+06 

-1.28E+06 

120 

36 

2.213 

9.429 

1.43E+06 

-1.43E+06 

120 

36 

2.213 

10.29 

i;56E+06 

-1.55E+06 

120 

36 

2.213 

11.14 

1.64E+06 

-1.63E+06 

120 

36 

2.213 

12 

1.69E+06 

-1.68E+06 

120 

36 

2.213 

12.86 

1.69E+06 

-1.68E+06 

120 

42 

2.049 

8.571 

1.37E+06 

-1.37E+06 

120 

42 

2.049 

9.429 

1.53E+06 

-1.53E+06 

120 

42 

2.049 

10.29 

1.67E+06 

-1.65E+06 

120 

42 

2.049 

11.14 

1.76E+06 

-1.74E+06 

120 

42 

2.049 

12 

1.81  E+06 

-1.79E+06 

120 

42 

2.049 

12.86 

1.81  E+06 

-1.79E+06 

120 

48 

1.917 

8.571 

1.45E+06 

-1 .45E+06 

120 

48 

1.917 

9.429 

1.62E+06 

-1.61  E+06 

120 

48 

1.917 

10.29 

1 .76E+06 

-1.75E+06 

120 

48 

1.917 

11.14 

1.86E+06 

-1.84E+06 

120 

48 

1.917 

12 

1.91  E+06 

-1.89E+06 

120 

48 

1.917 

12.86 

1.91  E+06 

-1.90E+06 

120 

54 

1.807 

8.571 

1.53E+06 

-1.52E+06 

120 

54 

1.807 

9.429 

1.71  E+06 

-1 .69  E+06 

120 

54 

1.807 

10.29 

1.85E+06 

-1 .83  E+06 

120 

54 

1.807 

11.14 

1.96E+06 

-1.92E+06 

120 

54 

1.807 

12 

2.01  E+06 

-1.97E+06 

120 

54 

1.807 

12.86 

2.01  E+06 

-1.97E+06 

120 

60 

1.714 

8.571 

1.59  E+06 

-1.58E+06 

120 

60 

1.714 

9.429 

1 .79E+06 

-1 .77E+06 

120 

60 

1.714 

10.29 

1 .94E+06 

-1.91  E+06 

120 

60 

1.714 

11.14 

2.05E+06 

-1.99E+06 

120 

60 

1.714 

12 

2.10E+06 

-2.03E+06 

120 

60 

1.714 

12.86 

2.10E+06 

-2.03E+06 

140 

35 

2.828 

0 

-2.66E+06 

2.68E+06 

1.87E+07 

-1.87E+07 

140 

35 

2.828 

-10.5 

-2.51  E+06 

2.52E+06 

1.66E+07 

-1.59E+07 

140 

35 

2.828 

-21 

-2.13E+06 

2.14E+06 

1.68E+07 

-1 .60E+07 

140 

35 

2.828 

-31.5 

-1.70E+06 

1.70E+06 

1 .70E+07 

-1.60E+07 

140 

35 

2.828 

-42 

-1.39E+06 

1.39E+06 

1.54E+07 

-1.79E+07 

140 

35 

2.828 

-52.5 

-1.30E+06 

1.32E+06 

1.63E+07 

-1.89E+07 

140 

42 

2.582 

0 

-2.65E+06 

2.68E+06 

2.68E+07 

-2.68E+07 

140 

42 

2.582 

-9.8 

-2.52E+06 

2.54E+06 

2.46E+07 

-2.60E+07 

140 

42 

2.582 

-19.6 

-2.15E+06 

2.16E+06 

2.49E+07 

-2.06E+07 

140 

42 

2.582 

-29.4 

-1.71  E+06 

1.72E+06 

2.52E+07 

-2.07E+07 

140 

42 

2.582 

-39.2 

-1.37E+06 

1 .37 E+06 

2.16E+07 

-2. 66  E+07 

140 

42 

2.582 

-49 

-1 .22E+06 

1.18E+06 

2.26E+07 

-2. 72  E+07 

140 

49 

2.39 

0 

-3.00E+06 

2.61  E+06 

3.51  E+07 

-3.51  E+07 

140 

49 

2.39 

-9.1 

-2.82E+06 

2.48E+06 

3.35E+07 

-3.41  E+07 

140 

49 

2.39 

-18.2 

-2.11  E+06 

2. 12  E+06 

3.39E+07 

-3.14E+07 

140 

49 

2.39 

-27.3 

-1.68E+06 

1.69E+06 

3.44E+07 

-2.35E+07 

140 

49 

2.39 

-36.4 

-1.31  E+06 

1.31  E+06 

3.46E+07 

-3.61  E+07 

140 

49 

2.39 

-45.5 

-1.10E+06 

1.05E+06 

2.83E+07 

-3.71  E+07 

140 

56 

2.236 

0 

-2.93E+06 

3.00E+06 

4.28E+07 

-4.28E+07 

140 

56 

2.236 

00 

1 

-2.76E+06 

2.35E+06 

4.27E+07 

-4.15E+07 

140 

56 

2.236 

-16.8 

-2.32E+06 

2.03E+06 

4.33E+07 

-3.83E+07 

140 

56 

2.236 

-25.2 

-1.74E+06 

1.61  E+06 

4.39E+07 

-4.27E+07 

140 

56 

2.236 

-33.6 

-1.21  E+06 

1.21  E+06 

4.43E+07 

-4.59E+07 

140 

56 

2.236 

-42 

-9.65E+05 

9.32E+05 

4.43E+07 

-4.73E+07 

140 

63 

2.108 

0 

-2.85E+06 

2.88E+06 

4.89E+07 

-4.89E+07 

140 

63 

2.108 

-7.7 

-2.65E+06 

2.73E+06 

5.16E+07 

-4.75E+07 

140 

63 

2.108 

-15.4 

-2.26E+06 

2.30E+06 

5.23E+07 

-4.38E+07 

140 

63 

2.108 

-23.1 

-1 .73E+06 

1.72E+06 

5.31  E+07 

-5.13E+07 

140 

63 

2.108 

-30.8 

-1 .20E+06 

1 .20E+06 

5.37E+07 

-5.51  E+07 

140 

63 

2.108 

-38.5 

-8.32E+05 

8.35E+05 

5.39E+07 

-6.06E+07 

140 

70 

2 

0 

-2.71  E+06 

2.72E+06 

5.88E+07 

-5.88E+07 

140 

70 

2 

-7 

-2.59E+06 

2.58E+06 

5.94E+07 

-5.14E+07 

140 

,  70 

2 

-14 

-2.16E+06 

2.22E+06 

6.03E+07 

-5.85E+07 

140 

70 

2 

-21 

-1 .68E+06. 

1 .71  E+06 

6.1 3  E+07 

-5.90E+07 

140 

70 

2 

-28 

-1 .28E+05 

1.19E+06 

6.21  E+07 

-6.33E+07 

140 

70 

2 

-35 

-1.73E+05 

7.69E+05 

6.26E+07 

-6.97E+07 

91 
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imm 
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F_tr_sag  (N) 

F_tr_ho.q  (N 

140 

35 

2.828 

10 

1.83E+06 

-1.83E+06 

140 

35 

2.828 

11 

1.97E+06 

-1.97E+06 

140 

35 

2.828 

12 

2.04E+06 

-2.04E+06 

140 

35 

2.828 

13 

2.11E+06 

-2.04E+06 

140 

35 

2.828 

14 

2.31  E+06 

-1.97E+06 

140 

35 

2.828 

15 

2.49  E+06 

-1.83E+06 

140 

42 

2.582 

10 

1.99E+06 

-1.98E+06 

140 

42 

2.582 

11 

2.14E+06 

-2.13E+06 

140 

42 

2.582 

12 

2.21  E+06 

-2.21  E+06 

140 

42 

2.582 

13 

2.29E+06 

-2.21  E+06 

140 

42 

2.582 

14 

2.50E+06 

-2.13E+06 

140 

42 

2.582 

15 

2.69E+06 

-1.98E+06 

140 

49 

2.39 

10 

2.12E+06 

-2.12E+06 

140 

49 

2.39 

11 

2.29E+06 

-2.28E+06 

140 

49 

2.39 

12 

2.37E+06 

-2.36E+06 

140 

49 

2.39 

13 

2.45E+06 

-2.36E+06 

140 

49 

2.39 

14 

2. 68 E+06 

-2.28E+06 

140 

49 

2.39 

15 

2.88E+06 

-2.1 2  E+06 

140 

56 

2.236 

10 

2.25E+06 

-2.24E+06 

140 

56 

2.236 

11 

2.42  E+06 

-2.41  E+06 

140 

56 

2.236 

12 

2.51  E+06 

-2.50E+06 

140 

56 

2.236 

13 

2.59E+06 

-2.49E+06 

140 

56 

2.236 

14 

2.84E+06 

-2.41  E+06 

140 

56 

2.236 

15 

3.05E+06 

-2.24E+06 

140 

63 

2.108 

10 

2.37E+06 

-2.36E+06 

140 

63 

2.108 

11 

2.55E+06 

-2.53E+06 

140 

63 

2.108 

12 

2.64E+06 

-2.62E+06 

140 

63 

2.108 

13 

2.72E+06 

-2.62E+06 

140 

63 

2.108 

14 

2.98E+06 

-2.53E+06 

140 

63 

2.108 

15 

3.21  E+06 

-2. 35  E+06 

140 

70 

2 

10 

2.48E+06 

-2.46E+06 

140 

70 

2 

11 

2.67E+06 

-2.64  E+06 

140 

70 

2 

12 

2.76E+06 

-2.74  E+06 

140 

70 

2 

13 

2.85E+06 

-2.74  E+06 

140 

70 

2 

14 

3.12E+06 

-2.64E+06 

140 

70 

2 

15 

3.35E+06 

-2.46E+06 

■SUM 

BIS 

160 

40 

3.232 

0 

-4.08E+06 

4.10E+06 

2.93E+07 

-2.93E+07 

160 

40 

3.232 

-12 

-3.81  E+06 

3.82E+06  . 

2.89E+07 

-2.82E+07 

160 

40 

3.232 

-24 

-3.12E+06 

2.83E+06 

2.77E+07 

-2.49E+07 

160 

40 

3.232 

-36 

-2.39E+06 

2.34E+06 

2.65E+07 

-2.46E+07 

160 

40 

3.232 

-48 

-1.97E+06 

1 .97E+06 

2.57E+07 

-2.88E+07 

160 

40 

3.232 

-60 

-2.03E+06 

2.03E+06 

2.75E+07 

-3.12E+07 

160 

48 

2.951 

0 

-4.1 3  E+06 

4.17E+06 

4.30E+07 

-4.30E+07 

160 

48 

2.951 

-11.2 

-3.87E+06 

3.90E+06 

4.24E+07 

-3.92E+07 

160 

48 

2.951 

-22.4 

-3.22E+06 

3.23E+06 

4.08E+07 

-3.62E+07 

160 

48 

2.951 

-33.6 

-2.46E+06 

2.46E+06 

3.89E+07 

-3.01  E+07 

160 

48 

2.951 

-44.8 

-1.94E+06 

1 .94E+06 

3.77E+07 

-4.26E+07 

160 

48 

2.951 

-56 

-1.86E+06 

1.86E+06 

3.76E+07 

-4.67E+07 

160 

56 

2.732 

0 

-4.06E+06 

4.12E+06 

5.76E+07 

-5.76E+07 

160 

56 

2.732 

-10.4 

-3.83E+06 

3.87E+06 

5.70E+07 

-4.94E+07 

160 

56 

2.732 

-20.8 

-3.22E+06 

3. 24 E+06 

5.49E+07 

-4.53E+07 

160 

56 

2.732 

-31.2 

■>2. 47  E+06 

2.47E+06 

5.65E+07 

-5.18E+07 

160 

56 

2.732 

-41.6 

-1.88E+06 

1 .88E+06 

5.05E+07 

-5.79E+07 

160 

56 

2.732 

-52 

-1.67E+06 

1 .67E+06 

4.99E+07 

-6.42E+07 

160 

64 

2.556 

0 

-4.27E+06 

3.98E+06 

7.22E+07 

-7.22E+07 

160 

64 

2.556 

-9.6 

-4.03E+06 

3.76E+06 

7.15E+07 

-7.10E+07 

160 

64 

2.556 

-19.2 

-3.14E+06 

3.17E+06 

6.91  E+07 

-6.83E+07 

160 

64 

2.556 

-28.8 

-2.43E+06 

2.44E+06 

7.25E+07 

-6.50E+07 

160 

64 

2.556 

-38.4 

-1 .80E+06 

1.80E+06 

7.62  E+07 

-7.36E+07 

160 

64 

2.556 

-48 

-1.46E+06 

1 .46  E+06 

6.18E+07 

-8.22E+07 

160 

72 

2.409 

0 

-4.11  E+06 

3.76E+06 

8.55E+07 

-8.55E+07 

160 

72 

2.409 

-8.8 

-3.90E+06 

3.56E+06 

8.47E+07 

-8.41  E+07 

160 

72 

2.409 

-17.6 

-3.32E+06 

3.04E+06 

8.21  E+07 

-8.09E+07 

160 

72 

2.409 

-26.4 

-2.33E+06 

2.35E+06 

8.82E+07 

-8.53E+07 

160 

72 

2.409 

-35.2 

-1.69E+06 

1.70E+06 

9.29E+07 

-8.85E+07 

160 

72 

2.409 

-44 

-1.26E+06 

1.10E+06 

9.63E+07 

-9. 93  E+07 

160 

80 

2.286 

0 

-3.87E+06 

3.47E+06 

9.63E+07 

-9.63E+07 

160 

80 

2.286 

-8 

-3.69E+06 

3.81  E+06 

9.55E+07 

-9.47E+07 

160 

80 

2.286 

-16 

-3.20E+06 

3.27E+06 

9.27E+07 

-9.12E+07 

160 

80 

2.286 

-24 

-2.48E+06 

2.51  E+06 

1 .02E+08 

-9.84E+07 

160 

80 

2.286 

-32 

-3.66E+05 

1.72E+06 

1 .08E+08 

-1.02E+08 

160 

80 

2.286 

-40 

-1.16E+05 

1.06E+06 

1.13E+08 

-1.14E+08 

nanai 

WSSEi 

F_tr_sag  (N) 

■zsm 

160 

40 

3.232 

11.43 

2.51  E+06 

-2.51  E+06 

160 

40 

3.232 

12.57 

2.58E+06 

-2.56E+06 

160 

40 

3.232 

13.71 

2.88E+06 

-2.49E+06 

160 

40 

3.232 

14.86 

3.1 5  E+06 

-3.15E+06 

160 

40 

3.232 

16 

3.37E+06 

-3.37E+06 

160 

40 

3.232 

17.14 

3.54E+06 

-3.53E+06 

160 

48 

2.951 

11.43 

2.73E+06 

-2.72E+06 

160 

48 

2.951 

12.57 

2.79E+06 

-2.78E+06 

160 

48 

2.951 

13.71 

3.12E+06 

-2.70E+06 

160 

48 

2.951 

14.86 

3.42E+06 

-2.50E+06 

160 

48 

2.951 

16 

3.66E+06 

-3.64E+06 

160 

48 

2.951 

17.14 

3.83E+06 

-3.82E+06 

160 

56 

2.732 

11.43 

2.92E+06 

-2.92E+06 

160 

56 

2.732 

12.57 

2.98E+06 

-2.98E+06 

160 

56 

2.732 

13.71 

3.34E+06 

-2.90E+06 

160 

56 

2.732 

14.86 

3.65E+06 

-2.68E+06 

160 

56 

2.732 

16 

3.91  E+06 

-3.89  E+06 

160 

56 

2.732 

17.14 

4.10E+06 

-4.08E+06 

160 

64 

2.556 

11.43 

3.09E+06 

-3.09E+06 

160 

64 

2.556 

12.57 

3.16E+06 

-3.16E+06 

160 

64 

2.556 

13.71 

3.54E+06 

-3.07E+06 

160 

64 

2.556 

14.86 

3.87E+06 

-2.84E+06 

160 

64 

2.556 

16 

4.14E+06 

-4.11  E+06 

160 

64 

2.556 

17.14 

4.34E+06 

-4.30E+06 

160 

72 

2.409 

11.43 

3.26E+06 

-3.25E+06 

160 

72 

2.409 

12.57 

3.32E+06 

-3.32  E+06 

160 

72 

2.409 

13.71 

3.72E+06 

-3.23E+06 

160 

72 

2.409 

14.86 

4.07E+06 

-2.98E+06 

160 

72 

2.409 

16 

4.35E+06 

-2.61  E+06 

160 

72 

2.409 

17.14 

4.56E+06 

-4.50E+06 

160 

80 

2.286 

11.43 

3.41  E+06 

-3. 39  E+06 

160 

80 

2.286 

12.57 

3.47E+06 

-3.46E+06 

160 

80 

2.286 

13.71 

3.89E+06 

-3.37E+06 

160 

80 

2.286 

14.86 

4.25E+06 

-3.12E+06 

160 

80 

2.286 

16 

4.55E+06 

-2.73E+06 

160 

80 

2.286 

17.14 

4.77E+06 

-4.69E+06 

180 

45 

3.637 

0 

-5.14E+06 

5.82E+06 

4.65E+07 

-4.65E+07 

180 

45 

3.637 

-13.5 

-4.86E+06 

4.87E+06 

4.21  E+07 

-4.04E+07 

180 

45 

3.637 

-27 

-4.15E+06 

4.15E+06 

4.24  E+07 

-3.93E+07 

180 

45 

3.637 

-40.5 

-3.31  E+06 

3.31  E+06 

4.26E+07 

-3.91  E+07 

180 

45 

3.637 

-54 

-2.70E+06 

2.70E+06 

3.83E+07 

-4.29E+07 

180 

45 

3.637 

-67.5 

-2.48E+06 

2.71  E+06 

4.03E+07 

-4.74E+07 

180 

54 

3.32 

0 

-5.12E+06 

5.96E+06 

6.69E+07 

-6.69E+07 

180 

54 

3.32 

-12.6 

-5.51  E+06 

4.88E+06 

6.24E+07 

-6.49  E+07 

180 

54 

3.32 

-25.2 

-4.18E+06 

4.19E+06 

6.28E+07 

-5.1 9  E+07 

180 

54 

3.32 

-37.8 

-3.34E+06 

3.34E+06 

6.31  E+07 

-5.13E+07 

180 

54 

3.32 

-50.4 

-2.66E+06 

2.66E+06 

6.31  E+07 

-6.33E+07 

180 

54 

3.32 

-63 

-2.33E+06 

2.45E+06 

5.60E+07 

-7.09  E+07 

180 

63 

3.073 

0 

-5.88E+06 

5.95E+06 

8.81  E+07 

-8.81  E+07 

180 

63 

3.073 

-11.7 

-5.52E+06 

4.75E+06 

8.49E+07 

-8.55E+07 

180 

63 

3.073 

-23.4 

-4.57E+06 

4.10E+06 

8.55E+07 

-7.92E+07 

180 

63 

3.073 

-35.1 

-3.40E+06 

3.27E+06 

8.60E+07 

-5.99E+07 

180 

63 

3.073 

-46.8 

-2.54E+06 

2.48E+06 

8.62E+07 

-8.60E+07 

180 

63 

3.073 

-58.5 

-2.11  E+06 

2.18E+06 

8.57E+07 

-9.72E+07 

180 

72 

2.875 

0 

-5.87E+06 

5.81  E+06 

1 .08E+08 

-1 .08E+08 

180 

72 

2.875 

-10.8 

-5.40E+06 

5.47E+06 

1.08E+08 

-1.05E+08 

180 

72 

2.875 

-21.6 

-4.54E+06 

3.90E+06 

1.09E+08 

-9.70E+07 

180 

72 

2.875 

-3.42E+06 

3.43E+06 

1.10E+08 

-1.08E+08 

180 

72 

2.875 

maam 

-2.43E+06 

2.43E+06 

1.10E+08 

-1.18E+08 

180 

72 

2.875 

-54 

-1.92E+06 

1 .39E+06 

1.10E+08 

-1.24E+08 

180 

81 

2.711 

0 

-5.70E+06 

5.56E+06 

1 .30E+08 

-1.30E+08 

180 

81 

2.711 

-9.9 

-5.40E+06 

5.26E+06 

1.31  E+08 

-1 .21  E+08 

180 

81 

2.711 

-19.8 

-4.40E+06 

4.46E+06 

1 .32E+08 

-1.12E+08 

180 

81 

2.711 

-29.7 

-3.37E+06 

3.39E+06 

1.33E+08 

-1.30E+08 

180 

81 

2.711 

-39.6 

-2.36E+06 

2.37E+06 

1.33E+08 

-1 .42  E+08 

180 

81 

2.711 

-49.5 

-1.69E+06 

1.69  E+06 

1 .33E+08 

-1.50E+08 

180 

90 

2.571 

0 

-5.21  E+06 

5.22E+06 

1 .49  E+08 

-1 .49  E+08 

180 

90 

2.571 

-9 

-5.16E+06 

4.96E+06 

1 .50E+08 

-1.31  E+08 

180 

90 

2.571 

-18 

-4.46E+06 

4.26E+06 

1.52E+08 

-1.48E+08 

180 

90 

2.571 

-27 

-3.25E+06 

3.29E+06 

1.53E+08 

-1 .49  E+08 

180 

90 

2.571 

-36 

-2.28E+06 

2.30E+06 

1.54E+08 

-1.64E+08 

180 

90 

2.571 

-45 

-3.27E+05 

1 .52E+06 

1.54E+08 

-1.73E+08 
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mm 

180 

45 

3.637 

12.86 

3.27E+06 

-3.08E+06 

180 

45 

3.637 

14.14 

3.69E+06 

-3.69E+06 

180 

45 

3.637 

15.43 

4.05E+06 

-4.04E+06 

180 

45 

16.71 

4.32E+06 

-4.31  E+06 

180 

45 

18 

4.49E+06 

-4.49E+06 

180 

45 

3.637 

19.29 

4.48E+06 

-4.56E+06 

180 

54 

3.32 

12.86 

3.55E+06 

-3.35E+06 

180 

54 

3.32 

14.14 

4.00E+06 

-4.00E+06 

180 

54 

3.32 

15.43 

4.39E+06 

-4.38E+06 

180 

54 

3.32 

16.71 

4.68E+06 

-4.67  E+06 

180 

54 

3.32 

18 

4.87E+Q6 

-4.86E+06 

180 

54 

3.32 

19.29 

4.85E+06 

-4.93E+06 

180 

63 

3.073 

12.86 

3.79E+06 

-3.59E+06 

180 

63 

3.073 

14.14 

4.28E+06 

-3.41  E+06 

180 

63 

3.073 

15.43 

4.69E+06 

-4.68E+06 

180 

63 

3.073 

5.00E+06 

-4.99E+06 

180 

63 

3.073 

18 

5.21  E+06 

-5.19E+06 

180 

63 

3.073 

19.29 

5.19E+06 

-5.27E+06 

180 

72 

2.875 

12.86 

4.02E+06 

-3.80E+06 

180 

72 

2.875 

14.14 

4.53E+06 

-3.61  E+06 

180 

72 

2.875 

15.43 

4.97E+06 

-4.96E+06 

180 

72 

2.875 

16.71 

5.30E+06 

-5.28E+06 

180 

72 

2.875 

18 

5.51  E+06 

-5.49E+06 

180 

72 

2.875 

19.29 

5.49E+06 

-5.58E+06 

180 

81 

2.711 

12.86 

4.22E+06 

-4.00E+06 

180 

81 

2.711 

14.14 

4.77E+06 

-3.80E+06 

180 

81 

2.711 

15.43 

-5.21  E+06 

180 

81 

2.711 

16.71 

5.57E+06 

-5.55E+06 

180 

81 

2.711 

18 

5.80E+06 

-5.77E+06 

180 

81 

2.711 

19.29 

5.78E+06 

-5.86E+06 

180 

90 

2.571 

12.86 

4.41  E+06 

-4.1 9  E+06 

180 

90 

2.571 

14.14 

4.98E+06 

-3.98E+06 

180 

90 

2.571 

15.43 

5.46E+06 

-5.44E+06 

180 

90 

2.571 

16.71 

5.83E+06 

-5.79E+06 

180 

90 

2.571 

18 

6.06E+06 

-6.02E+06 

180 

90 

2.571 

19.29 

6.04E+06 

-6.10E+06 

wmvm 

Hi 

200 

50 

4.041 

0 

-7.41  E+06 

7.44E+06 

6.70E+07 

-6.70E+07 
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22.29 

8.92E+06 

-8.92E+06 

240 

60 

4.849 

24 

9.46E+06 

-9.46E+06 

240 

60 

25.71 

9.82E+06 

-9.82E+06 

240 

72 

4.426 

17.14 

7.83E+06 

-7.83E+06 

240 

72 

4.426 

18.86 

8.01  E+06 

-8.14E+06 

240 

72 

4.426 

20.57 

8.91  E+06 

-8.10E+06 

240 

72 

4.426 

22.29 

9.67E+06 

-7.71  E+06 

mm 

72 

4.426 

24 

1.03E+07 

-1.03E+07 

240 

72 

4.426 

25.71 

-1.06E+07 

240 

84 

4.098 

17.14 

8.38E+06 

-8.38E+06 

240 

84 

4.098 

18.86 

8.57E+06 

-8.72E+06 

240 

84 

4.098 

20.57 

9.53E+06 

-8.67E+06 

240 

84 

4.098 

22.29 

1.03E+07 

-8.25E+06 

240 

84 

4.098 

24 

1.10E+07 

-1.10E+07 

240 

84 

4.098 

25.71 

1.14E+07 

-1.14E+07 

240 

96 

3.833 

17.14 

8.89E+06 

-8. 88  E+06 

240 

96 

3.833 

18.86 

9.07E+06 

-9.24  E+06 

240 

96 

3.833 

20.57 

1.01E+07 

-9.19E+06 

240 

96 

3.833 

22.29 

1.10E+07 

-8.74E+06 

240 

96 

3.833 

24 

1.16E+07 

-7.92E+06 

240 

96 

3.833 

25.71 

1.21E+07 

-1.20E+07 

240 

108 

3.614 

17.14 

9.35E+06 

-9.34E+06 

240 

108 

3.614 

18.86 

9.54E+06 

-9.72E+06 

240 

108 

3.614 

20.57 

1.06E+07 

-9.67E+06 

240 

108 

3.614 

22.29 

1.15E+07 

-9.20E+06 

240 

108 

3.614 

24 

1.22E+07 

-8.34E+06 

240 

108 

3.614 

25.71 

1 .27E+07 

-1.27E+07 

240 

120 

3.429 

17.14 

9.78E+06 

-9.77E+06 

240 

120 

3.429 

18.86 

9.97E+06 

-1.02E+07 

240 

120 

3.429 

20.57 

1.11E+07 

-1.01E+07 

240 

120 

3.429 

22.29 

1 .20E+07 

-9.62E+06 

240 

120 

3.429 

24 

1.28E+07 

-8.72E+06 

240 

120 

3.429 

25.71 

1 .33E+07 

-1.32E+07 

260 

65 

5.253 

0 

-1 .43E+07 

1 .56E+07 

1.73E+08 

-1 .73E+08 

260 

65 

5.253 

-19.5 

-1 .35E+07 

1.35E+07 

1 .74E+08 

-1.67E+08 

260 

65 

5.253 

-39 

-1.13E+07 

1.13E+07 

1.65E+08 

-1.47E+08 

260 

65 

5.253 

-58.5 

-8.95E+06 

8.95E+06 

1 .67E+08 

-1.48E+08 

260 

65 

5.253 

00 

i 

-7.35E+06 

7.35E+06 

1.50E+08 

-1.76E+08 

260 

65 

5.253 

-97.5 

-7.06E+06 

7.29E+06 

1 .61 E+08 

-1 .83E+08 

260 

78 

4.795 

0 

-1.59E+07 

1 .60E+07 

2.60E+08 

-2.60E+08 

260 

78 

4.795 

-18.2 

-1.48E+07 

1.36E+07 

2.54E+08 

-2.33E+08 

260 

78 

4.795 

-36.4 

-1.21E+07 

1.15E+07 

2.43E+08 

-2.14E+08 

260 

78 

4.795 

-54.6 

-9.08E+06 

9.08E+06 

2.48E+08 

-1.88E+08 

260 

78 

4.795 

-72.8 

-7.23E+06 

7.23E+06 

2.51  E+08 

-2.47E+08 

260 

78 

4.795 

-91 

-6.57E+06 

6.59E+06 

2.21  E+08 

-2.74E+08 

260 

91 

4.439 

0 

-1.59E+07 

1 .60E+07 

3.47E+08 

-3.47E+08 

260 

91 

4.439 

-16.9 

-1.49E+07 

1 .33E+07 

3.40E+08 

-2.96E+08 

260 

91 

4.439 

-33.8 

-1.23E+07 

1.14E+07 

3.31  E+08 

-2.71  E+08 

260 

91 

4.439 

-50.7 

-9.13E+06 

8.95E+06 

3.39E+08 

-3.37E+08 

260 

91 

4.439 

-67.6 

-6.94E+06 

6.94E+06 

3.44E+08 

-3.36E+08 

260 

91 

4.439 

-84.5 

-5.93E+06 

5.86E+06 

2.88E+08 

-3.76E+08 

260 

104 

4.153 

0 

-1 .47E+07 

1 .56E+07 

4.32E+08 

-4.32E+08 

260 

104 

4.153 

-15.6 

-1 .48E+07 

1.47E+07 

4.16E+08 

-4.24E+08 

260 

104 

4.153 

-31.2 

-1 .22E+07 

1.09E+07 

4.24E+08 

-4.02E+08 

260 

104 

4.153 

-46.8 

-9.18E+06 

8.61  E+06 

4.33E+08 

-4.28E+08 

260 

104 

4.153 

-62.4 

-6.53E+06 

6.53E+06 

4.41  E+08 

-4.60E+08 

260 

104 

4.153 

-78 

-5.16E+06 

5.16E+06 

4.44E+08 

-4.83E+08 

260 

117 

3.915 

0 

-1 .44E+07 

1.49E+07 

5.08E+08 

-5.08E+08 

260 

117 

3.915 

mu 

-1 .44E+07 

1.41  E+07 

5.04E+08 

-4.98E+08 

260 

117 

3.915 

-28.6 

-1 .22E+07 

1.19E+07 

5.13E+08 

-4.74E+08 

260 

117 

3.915 

-42.9 

-9.06E+06 

9.08E+06 

5.25E+08 

-5.17E+08 

260 

117 

3.915 

-57.2 

-6.35E+06 

6.36E+06 

5.36E+08 

-5.57E+08 

260 

117 

3.915 

-71.5 

-4.55E+06 

4.55E+06 

5.42E+08 

-5.86E+08 

260 

130 

3.714 

0 

-1.38E+07 

1 .40E+07 

5.67E+08 

-5.67E+08 

260 

130 

3.714 

-13 

-1.18E+07 

1 .33E+07 

5.82E+08 

-5.57E+08 

260 

130 

3.714 

-26 

-1 .02E+07 

1.19E+07 

5.93E+08 

-5.30E+08 

260 

130 

3.714 

-39 

-9.12E+06 

9.16E+06 

6.07E+08 

-5.96E+08 

260 

130 

3.714 

-52 

-6.17E+06 

6.15E+06 

6.20E+08 

-6.44E+08 

260 

130 

3.714 

-65 

-7.55E+05 

4.07E+06 

6.30E+08 

KKi.l-JMJH 

mm 

EX53I 

BSOI 

F_tr_hofl  (N) 

260 

65 

5.253 

18.57 

8.34E+06 

260 

65 

5.253 

20.43 

9.41  E+06 

-9.41  E+06 

260 

65 

5.253 

22.29 

1.03E+07 

-1.03E+07 

260 

65 

5.253 

24.14 

1.10E+07 

-1.10E+07 

260 

65 

5.253 

26 

1.14E+07 

-1.14E+07 

260 

65 

5.253 

27.86 

1.19E+07 

-1.16E+07 

260 

78 

4.795 

18.57 

9.04E+06 

-9.25E+06 

260 

78 

4.795 

20.43 

1.02E+07 

-9.26E+06 

260 

78 

4.795 

22.29 

1.12E+Q7 

-1.12E+07 

260 

78 

4.795 

24.14 

1.19E+07 

-1.19E+07 

260 

78 

4.795 

26 

1.24E+07 

-1 .24E+07 

260 

78 

4.795 

27.86 

1.29E+07 

-1.26E+07 

260 

91 

4.439 

18.57 

9.67E+06 

-9.91  E+06 

260 

91 

4.439 

20.43 

1 .09E+07 

-9.91  E+06 

260 

91 

4.439 

22.29 

1.20E+07 

-1.19E+07 

260 

91 

4.439 

24.14 

1.27E+07 

-1.27E+07 

260 

91 

4.439 

26 

1.32E+07 

-1.32E+07 

260 

91 

4.439 

27.86 

1.38E+07 

-1.34E+07 

260 

104 

4.153 

18.57 

1.02E+07 

-1 .05E+07 

260 

104 

4.153 

20.43 

1.16E+07 

-1.05E+07 

260 

104 

4.153 

22.29 

1.27E+07 

-9.96E+06 

260 

104 

4.153 

24.14 

1.35E+07 

-1.35E+07 

260 

104 

4.153 

26 

1.40E+07 

-1.40E+07 

260 

104 

4.153 

27.86 

1.46E+07 

-1.42E+07 

260 

117 

3.915 

18.57 

1.08E+07 

-1.11E+07 

260 

117 

3.915 

20.43 

1.22E+07 

-1.11E+07 

260 

117 

3.915 

22.29 

1.33E+07 

-1.05E+07 

260 

117 

3.915 

24.14 

1.42E+07 

-1.42E+07 

260 

117 

3.915 

26 

1.47E+07 

-1.47E+07 

260 

117 

3.915 

27.86 

1.53E+07 

-1.49E+07 

260 

130 

3.714 

18.57 

1.13E+07 

-1.16E+07 

260 

130 

3.714 

20.43 

1 .27E+07 

-1.16E+07 

260 

130 

3.714 

22.29 

1 .39E+07 

-1.10E+07 

260 

130 

3.714 

24.14 

1.48E+07 

-1.48E+07 

260 

130 

3.714 

26 

1 .54E+07 

-1.54E+07 

260 

130 

3.714 

27.86 

1.60E+07 

-1.56E+07 

Bs 


Xs 


280 

70 

5.657 

0 

-1.86E+07 

1 .86E+07 

2.37E+08 

-2. 37  E+08 

280 

70 

5.657 

-21 

-1 .73E+07 

1 .73E+07 

2.21  E+08 

-2.10E+08 

280 

70 

5.657 

-42 

-1 .41 E+07 

1.34E+07 

2.18E+08 

-1 .88E+08 

280 

70 

5.657 

-63 

-1.08E+07 

1.08E+07 

2.15E+08 

-1 .86E+08 

280 

70 

5.657 

-84 

-8.67E+06 

8.78E+06 

2.00E+08 

-2.26E+08 

280 

70 

5.657 

-105 

-9.17E+06 

9.17E+06 

2.05E+08 

-2.42E+08 

280 

84 

5.164 

0 

-1.89E+07 

1.89E+07 

3.41  E+08 

-3.41  E+08 

280 

84 

5.164 

-19.6 

-1.77E+07 

1 .77E+07 

3.26E+08 

-3.31  E+08 

280 

84 

5.164 

-39.2 

-1.46E+07 

1.34E+07 

3.23E+08 

-2.69E+08 

280 

84 

5.164 

-58.8 

-1.11  E+07 

1 .08 E+07 

3.18E+08 

-2.30E+08 

280 

84 

5.164 

-78.4 

-8.67E+06 

8.67E+06 

3.14E+08 

-3.35E+08 

280 

84 

5.164 

-98 

-8.39E+06 

8.39E+06 

2.85E+08 

-3.63E+08 

280 

98 

4.781 

0 

-1.86E+07 

1.88E+07 

4.50E+08 

-4.50E+08 

280 

98 

4.781 

-18.2 

-1.82E+07 

1 .76E+07 

4.42E+08 

-4.37E+08 

280 

98 

4.781 

-36.4 

-1.47E+07 

1 .47E+07 

4.39  E+08 

-4.05E+08 

280 

98 

4.781 

-54.6 

-1.12E+07 

1.12E+07 

4.33E+08 

-4.42E+08 

280 

98 

4.781 

-72.8 

-8.43E+06 

8.14E+06 

4.27E+08 

-4.58E+08 

280 

98 

4.781 

-91 

-7.48E+06 

7.48E+06 

4.23E+08 

-5.01  E+08 

280 

112 

4.472 

0 

-1.84E+07 

1.82E+07 

5.51  E+08 

-5.51  E+08 

280 

112 

4.472 

-16.8 

-1.80E+07 

1.71  E+07 

5.62E+08 

-5.36E+08 

280 

112 

4.472 

-33.6 

-1.50E+07 

1.44E+07 

5.57E+08 

-4.97E+08 

280 

112 

4.472 

-50.4 

-1.11  E+07 

1.10E+07 

5.71  E+08 

-5.47E+08 

280 

112 

4.472 

-67.2 

-8.09E+06 

8.09E+06 

6.05E+08 

-5.86E+08 

280 

112 

4.472 

00 

1 

-6.56E+06 

6.56E+06 

5.37E+08 

-6.46E+08 

280 

126 

0 

-1.79E+07 

1.86E+07 

6.75E+08 

-6.75E+08 

280 

126 

4.216 

-15.4 

-1.70E+07 

1.76E+07 

6.75E+08 

-6.18E+08 

280 

126 

4.216 

-30.8 

-1 .48E+07 

1 .39 E+07 

6.70E+08 

-5.73E+08 

280 

126 

4.216 

-46.2 

-1.11  E+07 

1.07E+07 

6.95E+08 

-6.57E+08 

280 

126 

4.216 

-61.6 

•  -7.56E+06 

7.67E+06 

7.39E+08 

-7.10E+08 

280 

126 

4.216 

-77 

-5.68E+06 

5.68E+06 

7.73E+08 

-7.86E+08 

280 

140 

4 

0 

-1.71  E+07 

1 .76E+07 

7.73E+08 

-7.73E+08 

280 

140 

4 

-14 

-1.49E+07 

1 .67E+07 

7.72E+08 

-7.69E+08 

280 

140 

4 

-28 

-1.29  E+07 

1 .43E+07 

7.68E+08 

-7.61  E+08 

280 

140 

4 

-42 

-1.08E+07 

1.10E+07 

8.07E+08 

-7.93E+08 

280 

140 

4 

-56 

-7.54E+06 

7.56E+06 

8.60E+08 

-8.21  E+08 

280 

140 

4 

-70 

-4.74E+06 

4.76E+06 

9.05E+08 

-9.12E+08 

Iff?! 

ml 

| 

1 

1 1 

280 

70 

5.657 

20 

1.05E+07 

-1.05E+07 

280 

70 

5.657 

22 

1.16E+07 

-1.16E+07 

280 

70 

5.657 

24 

1 .25E+07 

-1.25E+07 

280 

70 

5.657 

26 

1.30E+07 

-1.30E+07 

280 

70 

5.657 

28 

1.37E+07 

-1.32E+07 

280 

70 

5.657 

30 

1.47E+07 

-1.31  E+07 

280 

84 

5.164 

20 

1.13E+07 

-1.13E+07 

280 

84 

5.164 

22 

1.26E+07 

-1.26E+07 

280 

84 

5,164 

24 

1.35E+07 

-1.35E+07 

280 

84 

5,164 

26 

1.41  E+07 

-1.41  E+07 

280 

84 

5.164 

28 

1.49E+07 

-1.43E+07 

280 

84 

5.164 

30 

1.59  E+07 

-1.42E+07 

280 

98 

4.781 

20 

1.21  E+07 

-1.12E+07 

280 

98 

4.781 

22 

1 .35  E+07 

-1.35E+07 

280 

98 

4.781 

24 

1.45E+07 

-1.45E+07 

280 

98 

4.781 

26 

1.51  E+07 

-1.51  E+07 

280 

98 

4.781 

28 

1 .59E+07 

-1.53E+07 

280 

98 

4.781 

30 

1 .70E+07 

-1 .52  E+07 

280 

112 

4.472 

20 

1.28E+07 

-1.19E+07 

280 

112 

4.472 

22 

1 .43  E+07 

-1.43E+07 

280 

112 

4.472 

24 

1 .53  E+07 

-1.53E+07 

280 

112 

4.472 

26 

1.60E+07 

-1.60E+07 

280 

112 

4.472 

28 

1 .68  E+07 

-1.63E+07 

280 

112 

4.472 

30 

1.80E+07 

-1.61  E+07 

280 

126 

4.216 

20 

1 .35  E+07 

-1.25E+07 

280 

126 

4.216 

22 

1.50E+07 

-1.19E+07 

280 

126 

4.216 

24 

1.61  E+07 

-1.61  E+07 

280 

126 

4.216 

26 

1.68  E+07 

-1.68E+07 

280 

126 

4.216 

28 

1.77E+07 

-1.71  E+07 

280 

126 

4.216 

30 

1.90E+07 

-1.69E+07 

280 

140 

4 

20 

1.41  E+07 

-1.31  E+07 

280 

140 

4 

22 

1 .57E+07 

-1.25E+07 

280 

140 

4 

24 

1.69  E+07 

-1.69E+07 

280 

140 

4 

26 

1.76E+07 

-1.76E+07 

280 

140 

4 

28 

1.85E+07 

-1.79E+07 

280 

140 

4 

30 

1.98  E+07 

-1.76E+07 

Pi 

Xs  (m) 

'  1  .  ;;1|:  JlfTB 

300 

75 

6.061 

0 

-2.15E+07 

2.27  E+07 

3.02E+08 

-3.02E+08 

300 

75 

6.061 

-22.5 

-2.10E+07 

2.02E+07 

2.93E+08 

-2.91  E+08 

300 

75 

6.061 

-45 

-1.68E+07 

1.68E+07 

2.82E+08 

-2.52E+08 

300 

75 

6.061 

-67.5 

-1.31E+07 

1.31  E+07 

2.83E+08 

-2.52E+08 

300 

75 

6.061 

-90 

-1.08E+07 

1.08E+07 

2.56E+08 

-2.95E+08 

300 

75 

6.061 

-112.5 

-1.07E+07 

1 .07E+07 

2.69E+08 

-3.02E+08 

300 

90 

5.533 

0 

-2.32E+07 

2. 32  E+07 

4.30E+08 

-4.30E+08 

300 

90 

5.533 

-21 

-2.16E+07 

2.04E+07 

4.30E+08 

-4.15E+08 

300 

90 

5.533 

-42 

-1.71E+07 

1.71  E+07 

4.14E+08 

-3.32E+08 

300 

90 

5.533 

-63 

-1.33E+07 

1.33E+07 

4.20E+08 

-3.30E+08 

300 

90 

5.533 

-84 

-1.06E+07 

1 .06E+07 

3.85E+08 

-4.39E+08 

300 

90 

5.533 

-105 

-9.90E+06 

9.90E+06 

3.84E+08 

-4.50E+08 

300 

105 

5.122 

0 

-2.27E+07 

2.32E+07 

5.86E+08 

-5.86E+08 

300 

105 

5.122 

-19.5 

-2.17E+07 

2.01  E+07 

5.79E+08 

-5.38E+08 

300 

105 

5.122 

-39 

-1.79E+07 

1 .70E+07 

5.55E+08 

KEE33 

300 

105 

5.122 

-58.5 

-1.32E+07 

1 .33E+07 

5.75E+08 

-5.72E+08 

300 

105 

5.122 

00 

r". 

■ 

-1.02E+07 

1 .02E+07 

5.91  E+08 

-6.01  E+08 

300 

105 

5.122 

-97.5 

-8.90  E+06 

8.90E+06 
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Appendix  C.  Wave  Height  Matching  Reliability  Index  of  Five 


Where  the  required  variables  in  this  algorithm  are  contained  in  Appendix  A 
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Appendix  D.  Fit  Parameters  for  Analytical  Thesis  Model 


Actual  by  Predicted  Plot  Longitudinal 


Summary  of  Fit  Longitudinal  Troughing 
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Actual  by  Predicted  Plot  Longitudinal 


Summary  of  Fit  Longitudinal  Cresting 
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Actual  by  Predicted  Plot  Positive  Twisting 


Summary  of  Fit  Positive  Twisting 

RSquare  0.998925 
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Summary  of  Fit  Negative  Twisting 
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Summary  of  Fit  Transverse  Sagging 
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POSITIVE  NUMBER:  CONSTRAINT  SATISFIED.  |  THESE  VALUES  ARE  NORMALIZED 

NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  |  BETWEEN  +1.  AND  -1. 

1.000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 

-2.000  :  CONSTRAINT  SUPPRESSED. 


INITIAL  FRAME  ADEQUACY  PARAMETER  VALUES  -  MODULE  1  OF  SUBSTR. 
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INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  2  OF  SUBSTR. 


POSITIVE  NUMBER:  CONSTRAINT  SATISFIED.  |  THESE  VALUES  ARE  NORMALIZED 

NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  |  BETWEEN  +1.  AND  -1. 

1.000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 

-2.000  :  CONSTRAINT  SUPPRESSED. 

:  STRAKE  NOT  EVALUATED. 


INITIAL  FRAME  ADEQUACY  PARAMETER  VALUES  -  MODULE  2  OF  SUBSTR. 
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INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  1  OF  SUBSTR. 
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POSITIVE  NUMBER:  CONSTRAINT  SATISFIED.  I  THESE  VALUES  ARE  NORMALIZED 

NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  j  BETWEEN  +1.  AND  -1. 

1.000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 

-2.000  :  CONSTRAINT  SUPPRESSED. 

:  STRAKE  NOT  EVALUATED. 


INITIAL  FRAME  ADEQUACY  PARAMETER  VALUES  -  MODULE  2  OF  SUBSTR. 
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Analytical  Cresting 

INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  3  OF  SUBSTR . 


POSITIVE  NUMBER:  CONSTRAINT  SATISFIED.  |  THESE  VALUES  ARE  NORMALIZED 
NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  |  BETWEEN  +1.  AND  -1. 

1.000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 

-2.000  :  CONSTRAINT  SUPPRESSED. 

:  STRAKE  NOT  EVALUATED. 


INITIAL  FRAME  ADEQUACY  PARAMETER  VALUES  -  MODULE  3  OF  SUBSTR. 
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POSITIVE  NUMBER:  CONSTRAINT  SATISFIED.  |  THESE  VALUES  ARE  NORMALIZED 
NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  |  BETWEEN  +1.  AND  -1. 

1.000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 
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INITIAL  FRAME  ADEQUACY  PARAMETER  VALUES  -  MODULE  1  OF  SUBSTR. 
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POSITIVE  NUMBER:  CONSTRAINT  SATISFIED.  |  THESE  VALUES  ARE  NORMALIZED 

NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  |  BETWEEN  +1.  AND  -1. 

1.000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 

-2.000  :  CONSTRAINT  SUPPRESSED. 

—  :  STRAKE  NOT  EVALUATED. 
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POSITIVE  NUMBER!  CONSTRAINT  SATISFIED.  |  THESE  VALUES  ARE  NORMALIZED 

NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  |  BETWEEN  +1.  AND  -1. 

1.000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 

-2.000  :  CONSTRAINT  SUPPRESSED. 

:  STRAKE  NOT  EVALUATED. 
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INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  2  OF  SUBSTR. 


POSITIVE  NUMBER:  CONSTRAINT  SATISFIED.  |  THESE  VALUES  ARE  NORMALIZED 
NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  |  BETWEEN  +1.  AND  -1. 

1.000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 

-2.000  :  CONSTRAINT  SUPPRESSED. 

-.  STRAKE  NOT  EVALUATED. 
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Analytical  Transverse  Sagging 

INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  2  OF  SUBSTR.  2 

T_stRAKe”I_PCSF  J~PCCB  |  PCMY  |  PCSB  |  PYTF  |  PYTP  |  PYCF  |  PYCP  |  PSPBT  |  PSPBL  |  PFLB 


POSITIVE  NUMBER:  CONSTRAINT  SATISFIED.  1  THESE  VALUES  ARE  NORMALIZED 
NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  |  BETWEEN  +1.  AND  -1. 

1.000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 

-2.000  :  CONSTRAINT  SUPPRESSED. 

:  STRAKE  NOT  EVALUATED. 
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POSITIVE  NUMBER:  CONSTRAINT  SATISFIED.  |  THESE  VALUES  ARE  NORMALIZED 

NEGATIVE  NUMBER:  CONSTRAINT  VIOLATED.  |  BETWEEN  +1.  AND  -1. 

I. 000  :  CONSTRAINT  NOT  RELEVANT  OR  NULLIFIED  BY  USER. 

-2.000  ;  CONSTRAINT  SUPPRESSED. 

v  :  STRAKE  NOT  EVALUATED. 
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INITIAL  PANEL  ADEQUACY  PARAMETER  VALUES  -  MODULE  3  OF  SUBSTR. 
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